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EXECUTIVE  SUMMARY 


The  ‘funnel/gate’  system  is  a  developing  technology  for  passive  groundwater  plume 
management  and  treatment.  This  technology  uses  impermeable  ftmnel  walls  to  force  polluted 
groundwater  through  a  highly  permeable  zone  of  reactive  porous  media,  the  ‘gate,’  where 
contaminants  are  degraded  by  biotic  or  abiotic  heterogeneous  reactions.  The  purpose  of  this 
study  is  to  first  develop  and  demonstrate  analytical  flow  and  fate/transport  models  for 
funnel/gate  systems,  and  second,  to  develop  and  demonstrate  a  minimum  cost  design  model. 

A  series  of  equations  which  constitutes  an  analytical,  multicontaminant,  degradation 
and  transport  model  has  been  developed.  This  model  can  be  applied  in  iterative  manner  to 
determine  the  critical  hydraulic  residence  time,  T,  required  to  achieve  specified  concentrations 
of  contaminants  within  a  gate  (e.g.,  to  be  certain  that  the  water  quality  at  the  gate  exit  meets  a 
specified  water  quality  standard  i.e.,  1  pg/1  vinyl  chloride). 

From  an  existing  analytical  hydraulic  model  for  funnel/gate  systems  developed  by 
Christensen  and  Hatfield  (1994),  an  iterative  design  algorithm  was  developed  for  making  a 

funnel/gate  system  dimensional.  The  design  algorithm  requires  the  following  site  and  system 
data  to  execute: 

a)  the  critical  hydraulic  residence  time  in  the  gate; 

b)  the  transverse  width  of  groundwater  flow  to  be  intercepted  within  the  funnel,  B ; 

c)  the  saturated  thickness  of  the  phreatic  aquifer,  (j)i  at  the  entrance  of  the  funnel; 

d)  the  groundwater  flow  through  the  area  defined  by  width  B  and  flow  thickness,  (j), ; 

e)  the  aquifer  hydraulic  conductivity; 

f)  the  porosity  of  the  reactive  porous  media  within  the  gate;  and 

g)  the  hydraulic  conductivity  of  the  porous  media  within  the  gate. 

Results  of  two  FRAC-3D  numerical  validation  studies  are  also  presented.  Thus,  results 
of  the  numerical  validation  suggest  the  funneFgate  design  model  could  be  used  to  predimension 


a  system;  however,  subsequent  numerical  simulations  should  be  conducted  to  verify  the 
hydraulics  of  the  design. 

A  fiinnel/gate  cost  estimation  model  was  also  developed  to  identify  minimum  cost 
designs  subject  to  four  constraints:  one  that  defines  the  desired  hydraulic  retention  time  in  the 
gate,  one  that  defines  the  funnel  wall  length  in  terms  of  other  system  dimensions,  and  two  that 
serve  as  hydraulic  constraints.  Costs  considered  in  the  model  include  land  costs,  costs  per  unit 
length  of  funnel  walls,  costs  per  unit  length  of  gate,  costs  per  unit  width  of  gate,  and  costs 
associated  with  the  reactive  medium  used  in  the  gate.  As  a  result  of  recasting  into  a  Lagrangian 
optimization  model,  the  optimization  problem  was  reduced  to  that  of  solving  a  system  of  nine 

nonlinear  equations  with  nine  unknowns.  Recommendations  were  given  as  to  how  the  nine 
equations  could  be  solved. 

Finally  a  computer  program  was  developed  to  demonstrate  the  theory  that  was  developed. 
The  Funnel/Oate  Design  Model  (FGDM)  software  is  a  FORTRAN  program  developed  to 
ascertain  feasible  dimensions  of  funnel-and-gate  systems  and  also  identify  the  lowest  cost  design. 
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SECTION  I 


INTRODUCTION 


OBJECTIVE 

The  purpose  of  this  study  is  to  first  develop  and  demonstrate  analytical  flow  and 
fate/transport  models  for  funnel/gate  systems,  and  second,  to  develop  and  demonstrate  a 
minimum  cost  design  model. 

BACKGROUND 

The  ‘funnel/gate’  system  is  a  developing  technology  for  passive  groimdwater  plume 
management  and  treatment.  This  technology  uses  impermeable  funnel  walls  to  force  polluted 
groundwater  through  a  highly  permeable  zone  of  reactive  porous  media,  the**gate^^^where 
contaminants  are  degraded  by  biotic  or  abiotic  heterogeneous  reactions. 

Figure  1  illustrates  a  plan  view  of  a  funnel/gate  system  of  symmetrical  configuration. 
This  system  contains  matching  sets  of  funnel  walls.  The  first  set  of  funnel  walls  serves  to 
intercept  the  plume  and  direct  polluted  groundwater  into  the  gate;  whereas,  the  second  set 
functions  to  disperse  the  decontaminated  groundwater  as  it  exits  the  gate.  The  funnel  walls  may 
be  slurry  walls  or  sheet  piles  keyed  into  an  underlying  aquitard. 

As  suggested  above,  groundwater  decontamination  occurs  inside  the  gate.  Here, 
dissolved  contaminant  removal  is  achieved  through  heterogeneous  degradation  reactions  on  the 
gate  porous  media  alone  or  through  the  action  of  biotic/abiotic  processes  supported  in  the  gate 
(i.e.,  bioremediation  or  sparging). 

To  design  a  funnel/gate  system  it  is  essential  to  first  estimate  d  ’’critical  hydraulic 
residence  time  "for  groimdwater  flowing  through  the  gate.  The  duration  of  this  time  is  chosen  to 
achieve  "adequate  contact  time'^between  the  intercepted  groundwater  and  the  reactive  gate 
porous  media.  In  this  report,  adequate  contact  time  refers  to  the  minimum  time  needed  to  bring 
polluted  groundwater  in  contact  with  the  reactive  gate  porous  media  such  that  dissolved 
contaminants  are  reduced  to  desired  levels;  that  is,  to  be  certain  water  quality  at  the  gate  exit 
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Direction  of  Groundwater  Flow 
- >► 


Figure  1.  Conceptual  Plan  View  of  a  Funnel/Gate  System. 

meets  specified  water  quality  goals  (i.e.,  system  performance  standards  or  drinking  water 
standards  i.e.,  1  pg/1  vinyl  chloride). 

Several  pollutants  may  be  targeted  for  destruction  in  a  fimnel/gate.  For  each  target 
contaminant,  the  adequate  contact  time  is  determined.  For  each  contaminant,  this  is  achieved 
using  the  observed  maximum  plume  concentrations,  the  expected  degradation  rate  coefficient  in 
the  gate,  and  the  water  quality  goal  applied  to  water  leaving  the  gate.  For  a  plume  comprised  of 
several  target  contaminants  (i.e.,  PCE  and  TCE),  a  critical  hydraulic  residence  time  must  be 
determined.  Here,  the  critical  hydraulic  residence  time  equals  the  longest  adequate  contact  time 
in  the  gate  which  ensures  all  target  contaminants  meet  desired  water  quality  goals  at  the  gate  exit. 
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Once  the  critical  hydraulic  residence  time  in  gate  is  known,  the  funnel  and  gate  can  be 

dimensioned  following  the  approach  of  Christensen  and  Hatfield  (1994),  and  using  information 
on: 

a)  the  transverse  width  of  groundwater  flow  to  be  intercepted  within  the  funnel,  this 
is  usually  taken  to  be  the  width  of  the  plume; 

b)  the  saturated  thickness  of  the  phreatic  aquifer  at  the  entrance  of  the  funnel; 

c)  the  hydraulic  gradient  of  groundwater  flow; 

d)  the  aquifer  permeability; 

e)  the  porosity  of  the  reactive  porous  media  within  the  gate;  and 

f)  the  permeability  of  the  porous  media  in  the  gate. 

SCOPE 

This  report  describes  hydraulic  and  contaminant  fate  and  transport  theories  applicable  to 
designing  funnel/gate  systems.  In  addition,  the  report  describes  a  software  package, 

Funnel/Gate  Design  Model  (FGDM),  which  was  developed  to  identify  minimum  cost  design 
configurations.  In  Section  2,  contaminant  fate  and  transport  theories  are  presented  as  they  apply 
to  dissolved  pollutant  degradation  inside  a  reactive  gate.  Section  3  reviews  the  hydraulics  of 
funnel/gate  systems  as  previously  presented  by  Christensen  and  Hatfield  (1994).  In  Section  4,  a 
design  model  for  minimizing  fiinnel/gate  system  costs  is  developed.  Section  5  reviews  the 
FGDM  software  and  results  from  an  example  problem.  Finally,  Section  6  lists  the  FGDM 
FORTRAN  code  with  example  input  and  output  files. 
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SECTION  II 


CONTAMINANT  FATE  AND  TRANSPORT  THEORY 
-tlie  first  objective  of  this  chapter  is  to  present  contaminant  fate  and  transport  theories 
pertment  to  flmnel/gate  systems.  A  second  but  equally  important  objective  is  to  develop  an 
analytieal  contaminant  degradation  and  transport  model  for  the  gate. 

THEORY 

Contaminant  fate  and  transport  hasreceived  considerable  attention  in  gmundwater 
_  systems.  For  fimnel/gate  systems,  research  has  focused  on  die  fate  and  transport  of  groundwater 
conmnunants  mrgeted  for  destruction  inside  the  gate.  Much  of  this  research  has  examined 
red^tive  dehalogenation  of  chlorinated  allcenes  and  alhanes  using  zem  valent  iron  (Oillham  and 

^  Hannesm,  1994;  Schreier  and  Reinhard,  1994;  Wanen,  Arnold,  Bishop.  Lindholm,  and 

Betterton,  1995;  and  Bunts.  Camphell.  and  Manoranjan,  1995).  Gillham  and  O'Hannesin  (1994) 

and  Burns'  have  looked  at  the  degradation  of  tetrachloroethylene  (PCE).  trichloroethylene 

(TCE),  dichloroethylene  (DCE),  vinyl  ehloride  and  other  halogenated  aliphatics  using  zero 

valent  non.  A  fraction  of  this  degradation  oceurs  within  a  sequential  reductive  dehalogenation 

F^ess  that  transforms  PCE  to  TCE,  TCE  to  DCE,  and  DCE  to  vinyl  chloride  (Gillham  and 

O  Hannesm,  1995;  and  Burris')  .In  addition  to  this  sequence  of  reactions,  Gillham  and 

O’Hannesin  (1995)  and  Burris'  have  indicated  evidence  of  competing  processes  whereby  PCE 

TCE.  DCE  are  transformed  to  other  produets  such  as  PCE  degradation  to  dichloroacetylene  ’ 

TCE  to  chloroacetylene,  and  DCE  degradation  to  acetylene.  Figme  2  summarizes  these 

Parameters  A.,^,  etc.,  appearing  in  the  figure  correspond  to  degradation  rate 
coefficients. 


Bums,  D.R.,  Personal  Communication 
Drive,  Tyndall  AFB,  FL  32403,  1995. 


Armstrong  Laboratory,  AL/EQC,  139  Bames 
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I  Vinyl  CMoride 

Figure  2.  Conceptual  Model  of  Reactions  Involving  Tetrachloroethylene, 
Trichloroethylene,  Dichloroethylene,  and  Vinyl  Chloride. 
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To  develop  a  model  describing  the  contaminant  fate  and  transport  inside  the  gate  of  a 
funnel/gate  system,  it  is  necessary  to  adopt  a  conceptual  model  of  the  reaction  configuration. 
That  is,  a  conceptual  model  descriptive  of  the  reactions  encountered  with  a  dissolved 
contaminant  inside  the  gate.  Following  the  reaction  configuration  illustrated  in  Figure  2,  a 
conceptual  model  was  developed  to  generally  describe  not  only  the  sequential  degradation  of 
chlorinated  alkenes  but  also  other  contaminants  (i.e.,  nitrate).  Thus,  Figure  3  shows  contaminant 
C,  degrading  to  C2,  then  C2  to  C3 ,  then  C3  to  C4 ,  and  finally  the  degradation  of  C4.  If  it  is 
assumed  that  C„  C2,  C3,  and  C4  respectively  correspond  to  contaminants  PCE,  TCE,  DCE,  and 
vinyl  chloride,  then  the  reaction  conceptual  model,  as  illustrated  in  Figure  3,  is  sufficient  to 
describe  the  degradation  of  these  contaminants  in  the  gate. 

Inside  the  gate,  advection  (as  opposed  to  dispersion)  is  the  dominant  transport 
mechanism.  Using  Figure  3  and  assuming  the  gate  flow  field  to  be  one-dimensional,  governing 
steady-state  fate  and  transport  equations  can  be  written  in  terms  of  contaminant  concentration 
variables  C„  C2,  C3,  and  C4.  Thus, 


ac, 

V - 

dX 


(1) 


^2A^2 


(2) 
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Figxye  3. 


General  Conceptual  Model  of  the  Reaction  Configuration  Inside  the  Gate  of 
a  Funnel/Gate  System. 
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V 


(3) 


dCj 

W 


dC, 

V - 

dx 


^4A^4  ■'■  ^3A^3 


(4) 


where  v  =  unidirectional  pore  water  velocity  for  groundwater  flowing  through  the  gate;  X  =  the 
longitudinal  distance  measured  from  the  gate  entrance  that  runs  parallel  to  the  direction  of 
groundwater  flowthrough  the  gate  (see  Figure  1);  and  A, a,  Ajb,  A2A,  ^2b>  ^3a>  ^3b>  Ka’  ~  first 
order  degradation  coefficients  as  illustrated  in  Figure  3. 

Equations  (1)  through  (4)  are  simple  advective  transport  equations  which  assume 
unidirectional,  uniform,  and  constant  pore  velocities  throughout  the  gate.  These  equations 
assume  contaminant  concentrations  are  invariant  along  transverse  Y  and  Z  directions; 
consequently,  they  are  limited  to  describing  steady-state  concentration  profiles  of  C,,  C2,  C3,  and 
C4  in  the  gate  and  along  the  X  axis  alone. 

ANALYTICAL  MULTI-CONTAMINANT  DEGRADATION  AND  TRANSPORT  MODEL 

Analytical  solutions  to  Equations  (1)  through  (4)  can  be  obtained  under  the  assumption 
that  influent  concentrations  are  constant  and  known  for  C,,  C2,  C3,  and  C4  at  the  gate  entrance  (X 
=  0).  This  is  equivalent  to  assuming  boundary  conditions  that  specify  concentrations  variables 
C,,  C2,  C3,  and  C4  as  having  values  of  Cjo,  C20,  C30,  and  C40  respectively  at  the  X  =  0.  In  the  field, 
contaminant  concentrations  are  never  constant  at  the  gate  entrance,  but  vary  with  time  as  a  plume 
is  intercepted.  As  a  result,  it  will  be  necessary  in  practice  to  let  Cjo,  C20,  C30,  and  C40  equal 
maximum  plume  concentrations  of  four  contaminants  identified  in  a  plume  and  targeted  for 
destruction  inside  the  gate.  Under  this  approach,  a  critical  hydraulic  residence  is  found  which  is 
suitable  for  treating  maximum  observed  pollutant  concentrations. 
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Using  the  above  boundary  assumption,  that  Ci  (X  =  0)  =  C,o,  the  solution  to  Equation  (1) 
is; 

C,{T)  =  C,,-EXP[-jT]  (5) 

where 

j  ~  ^\A  ^  (^) 

and  T  =  the  hydraulic  residence  time  as  defined  by  Eq.  (7)  below: 

T=-  (7) 

V 

Equation  (5)  describes  the  steady-state  concentration  profile  of  contaminant  C,  throughout  the 
length  of  the  gate. 

Using  Equation  (5)  and  the  boundary  condition  defining  C2  (X  =  0)  =  C20,  the  solution  to 
Equation  (2)  is: 

KC 

C^iT)  =  C^o-EXPi- gT]  +- - ^±-[EXP[-jT]  -  EXP[- gT  ]  ]  (8) 

(g-j) 

where 

S  ~  ^2A  ^2B 

and 

K  =  (10) 

Equation  (8)  describes  the  steady-state  concentration  profile  of  contaminant  Cj  throughout  the 
length  of  the  gate  as  a  consequence  of  Cj  degradation  and  of  Cj  production  from  C,  degradation. 

Again,  Equation  (3)  is  solved  similarly,  using  Equation  (8)  and  the  above  stated  boundary 
condition  of  C3  (X  =  0)  =  C30.  Thus, 
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Where 


C^{T)  =  A-EXP[-  mT]  +  E-EXP[- jT]  +  D-EXP[-  gT] 


A  =  a 


f^20  ^  fKC,, 

is  -  j){^  -  j)  (g  -  j)(m  -  g) 


(g  -j){m  -  g) 


fKC,, 

(g  -  j){m  -  j) 


/^2o  _  fKC,, 

-  g)  (g-j)(fn-g) 


f=K 


Thus,  Equation  (11)  describes  the  steady-state  concentration  profile  of  contaminant  C3 
throughout  the  length  of  the  gate  as  both  a  consequence  of  C3  degradation  to  C4  and  of  C3 
production  from  Cj  degradation. 
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Finally,  substituting  Equation  (1 1)  into  Equation  (4)  and  using  the  above  stated  boundary 
condition  of  C4  (X  =  0)  =  C40,  the  following  function  is  obtained  for  C4  concentration 
distributions  inside  the  gate: 

P  A 

CAT)  =  G  EXP[  -  HT]  +  ^J--—EXP[-  mT] 

^  (H  -  m) 

+  ^?—EXP[-  jT]  +  ^^  -EXPi-  gT]  (17) 

{H  -  j)  {H  -  g) 


Equation  (17)  describes  steady-state  C4  concentration  profiles  that  reflect  degradation  and  C4 
production  from  C3  degradation. 

SUMMARY 

Equations  (5),  (8),  (1 1),  and  (17)  constitute  an  analytical,  multi-contaminant,  degradation 
and  transport  model.  This  model  can  be  applied  in  iterative  manner  to  determine  the  critical 
hydraulic  residence  time,  T,  required  to  achieve  specified  concentrations  of  contaminants  C,,  C2, 
C3,  and  C4  within  a  gate  (e.g.,  to  be  certain  that  the  water  quality  at  the  gate  exit  meets  a  specified 
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water  quality  standard  such  as  1  qg/l  vinyl  chloride).  To  use  this  transport  model,  the  following 
data  are  needed: 

a)  degradation  rate  coefficients  for  all  reactions  depicted  in  Figure  3; 

b)  specified  water  quality  goals  or  system  performance  standards  for  Ci,  C3,  and 
C4  (e.g.,  primary  drinking  water  standards  for  PCE,  TCE,  DCE,  and  vinyl 
chloride);  and 

c)  the  concentration  of  each  contaminant  at  the  entrance  of  the  gate,  C,o,  C20,  C30, 
and  C40.  Once  the  critical  hydraulic  residence  time  is  known,  the  hydraulic  model 
can  be  used  to  dimension  the  system. 
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SECTION  III 


FUNNEL/GATE  HYDRAULICS 

This  Section  seeks  to  fulfill  two  objectives.  The  first  is  to  review  elements  of  an  existing 
analytical  hydraulic  model  for  funnel/gate  systems  that  was  developed  by  Christensen  and 
Hatfield  (1994).  The  second  objective  is  to  compare  analytical  model  predictions  of  funnel/gate 
system  performance  to  those  given  by  a  complex  numerical  model,  FRAC-3D. 

THEORY 

Christensen  and  Hatfield  (1994)  presented  a  model  describing  the  steady-state  hydraulics 
of  a  funnel/gate  system.  Model  development  was  predicated  on  the  symmetrical  system 
configuration  illustrated  in  Figure  4.  This  figure  presents  a  conceptual  plan  view  of  a  funnel/gate 
system  about  to  capture  a  plume  of  width  B,,.  In  addition,  this  figure  shows  groundwater,  of  flow 
width  less-than-or-equal-to  ‘B,’  being  intercepted  from  an  undisturbed  flow  field  by  an  upstream 
funnel  of  entrance  width  B,  and  projected  wall  length,  L.  As  a  result  of  gravity,  groundwater  is 
forced  through  a  gate  of  width,  b,  and  length,  {.  The  hydraulic  residence  time  inside  the  gate  is 
T.  As  this  groundwater  exits  the  gate,  flow  diverges  in  a  downstream  funnel  identical  to  the 
upstream  funnel  and  is  returned  to  the  undisturbed  flow  field. 

Christensen  and  Hatfield  (1994)  characterized  groundwater  flow  up  gradient,  down 
gradient,  and  inside  the  funnel/gate  by  dividing  the  overall  flow  field  into  five  zones  (See  Figure 
5).  Zones  1  and  5  respectively  correspond  to  up  gradient  and  down  gradient  undisturbed 
groundwater  flow  fields.  Zone  2  is  a  converging  flow  zone  created  inside  the  upstream  funnefi 
whereas,  Zone  4  encompasses  a  region  of  divergent  flow  within  the  downstream  funnel.  Finally, 
there  is  Zone  3,  inside  the  gate,  where  flow  is  assumed  one-dimensional  and  horizontal. 

As  stated  above.  Zones  1  and  5  reflect  undisturbed  flow  fields;  however,  in  Zones  2, 3, 
and  4  there  are  changes  in  groundwater  potential  caused  by  a  fimnel/gate  system  that  must  be 
estimated.  Christensen  and  Hatfield  (1994)  developed  their  hydraulic  model  treating  the  aquifer 
as  an  isotropic  and  homogeneous  medium  underlain  by  a  horizontal  aquiclude.  In  addition,  they 
assumed  the  capillary  fnnge  was  minor,  and  that  groundwater  flow  was  undirectional,  horizontal, 
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Figure  4.  Plan  View  of  Symmetrical  Fuimel/Gate  System  Defining  Geometry  of  Flow 
Boundaries. 
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and  unconfined.  From  these  assumptions,  they  derived  the  following  functions  which 
characterize  groundwater  potentials  inside  zones  2, 3,  and  4,  respectively: 


<^2(^2) 


( 


2QL 

k(B  -  b) 


1 


b 

X2' 

B 

L] 

(21) 


k{B  -  b)  B 


2QX, 

k^b 


(22) 


+ 


2QL  ^^b  _  2^  _ 
k(B  -  b)  B  k^b 


2QL 

k{B  -  b) 


In 


/ 

1  - 

[1  --I 

\ 

B 

l] 

(23) 


where  (^>2,  <1)3,  and  =  groundwater  potentials  in  Zones  2,  3,  and  4,  respectively  (these  are  also 
equivalent  to  aquifer  saturated  thicknesses  in  these  zones  if  the  zero  datum  is  taken  to  be  the  top 
of  the  aquiclude);  Xj ,  X3 ,  and  X4  =  longitudinal  distances  from  the  upstream  ends  of  Zones  2,  3, 
and  4,  respectively  (see  Figure  4);  (})i  =  the  groundwater  potential  (or  saturated  aquifer  thickness) 
at  the  entrance  of  the  upstream  funnel  (where  X2  =  0  at  the  beginning  of  Zone  2),  k  —  the 
hydraulic  conductivity  of  aquifer;  k^  =  the  hydraulic  conductivity  of  the  gate;  and  Q  =  the  total 
rate  of  flow  passing  the  width  B  in  the  undisturbed  groundwater  flow  field. 

The  total  rate  of  flow  through  the  funnel/gate  system,  Q  is  estimated  using  the  aquifer 
permeability,  the  width  (  B)  of  groundwater  flow  intercepted,  the  potential  <})„  at  the  funnel 
entrance,  and  the  potential  gradient,  d(l)/dx,  of  the  undisturbed  flow  field  at  Xj  =  0.  Thus, 


Q  =  B^^k 


(24) 
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The  hydraulic  residence  time  for  the  groundwater  flowing  through  the  gate,  T,  is  given 
by: 


T  - 

Q 


(25) 


where  0  -  porosity  of  the  gate  and  (j)ga,e  =  the  average  groundwater  potential  in  the  gate,  which  is 
taken  here  to  equal  the  saturated  thickness  of  the  aquifer  inside  the  gate.  Design  parameters  Q 
and  6  are  usually  known  or  specified,  while  T  equals  the  critical  hydraulic  residence  time 
identified  by  the  transport  model.  Knowing  that  Q,  6,  and  T  are  constants,  it  follows  from 

Equation  (25)  that  the  gate  volume,  fb(t)g,,e,  is  also  constant  for  an  identified  critical  hydraulic 
residence  time. 

Generally,  it  may  be  assumed  that  the  width  of  intercepted  groundwater  flow  is  less-than- 
or-equal-to  B,  the  funnel  width,  and  that  B  ^  B  the  plume  width.  Under  this  assumption, 
Christensen  and  Hatfield  (1994)  found  that  the  projected  length  of  the  funnel  walls,  L  can  be 
estimated  from: 


k/k 

— ^  -  1 
g.  b/B 

2  Injb/B)  ^  j 

1  -  (b/B) 


(26) 


ITERATIVE  SOLUTION  ALGORITHM  TO  DIMENSION  A  FUNNEL/GATE  SYSTEM 

To  design  a  funnel/gate  system,  it  is  difficult  to  use  Equations  (25)  and  (26)  directly. 

This  is  because  gate  length  f  and  potential  are  generally  unknown,  and  also  because 
governing  Equations  (26),  (25),  and  (21)  or  (22)  must  be  satisfied  simultaneously.  As  a  result,  L, 
6,  and  <j)ga,e  are  estimated  using  an  iterative  approach. 
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The  design  process  begins  assuming  the  funnel  width  B  and  the  gate/funnel  width  ratio 
b/B  are  known;  therefore,  b  is  known.  In  practice,  a  suite  of  b/B  ratios  should  be  investigated. 
During  the  first  step  of  the  first  iteration,  the  gate  length  is  approximated  using  Equation  (25)  and 
substituting  (j)i  as  an  initial  approximation  for  (j)ga,e.  Thus, 


TQ 

06$, 


(27) 


The  next  step  uses  the  S  results  of  Equation  (27)  in  Equation  (26)  to  obtain  an  initial  estimate  of 
the  projected  funnel  wall  length,  L.  At  the  end  of  the  first  iteration,  initial  estimates  of  f  and  L 
have  been  made. 

During  the  first  step  of  the  second  iteration,  (Jigate  is  calculated  using  the  function  below; 


InA 

k(B  -  b)  B 


(28) 


Equation  (28)  can  be  derived  either  from  Equation  (21)  assuming  Xj  =  L  or  from  Equation  (22) 
assuming  Xj  =  0.  The  second  step  uses  the  most  recent  estimate  of  (})ga,e  in  Equation  (25)  to 
obtain  an  updated  value  for  the  gate  length,  H.  The  third  and  final  step  of  the  second  iteration 
uses  the  most  recent  estimate  of  H  in  Equation  (26)  to  obtain  a  new  calculation  of  the  projected 
funnel  wall  length,  L. 

Subsequent  iterations  repeat  steps  described  under  the  second  iteration;  however  during 
each  step  f ,  L,  and  (|)ga,e  estimates  generated  from  the  most  recent  iteration  are  used.  The  iterative 
process  continues  until  changes  in  {,  L,  and  (j)ga,e  values  are  less  than  1  percent  between  iterations. 

From  the  above  presentation,  it  may  be  concluded  that  the  hydraulic  model  can  be  used  to 
dimension  a  funnel/gate  system  given  the  following  information  is  available: 

a)  the  critical  hydraulic  residence  time  inside  the  gate; 

b)  the  transverse  width  of  groundwater  flow  to  be  intercepted  within  the  funnel,  B, 
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this  is  practically  taken  to  be  greater  than  the  width  of  the  plume; 

c)  the  saturated  thickness  of  the  phreatic  aquifer,  (j)]  at  the  entrance  of  the  funnel; 

d)  the  groundwater  flow  through  the  area  defined  by  width  B  and  flow  thickness,  (j)]; 

e)  the  aquifer  hydraulic  conductivity; 

f)  the  porosity  of  the  reactive  porous  media  within  the  gate;  and 

g)  the  hydraulic  conductivity  of  the  porous  media  within  the  gate. 

MODEL  VALIDATION 

To  validate  the  funnel/gate  hydraulic  model,  numerical  simulations  were  performed  to 
generate  flow  fields  for  several  fimnel/gate  configurations  dimensioned  using  the  above 
analytical  hydraulic  model.  The  numerical  simulations  were  conducted  using  FRAC-3D.  The 
velocity  fields  generated  from  the  numerical  simulations  were  then  used  to  create  figures 
illustrating  the  groundwater  flow  streamlines.  Figure  5  illustrates  conceptually  the  orientation  of 
a  fimnel/gate  system  in  the  numerical  model.  Shown  is  a  wide  funnel/gate  system.  The  gate 
width,  b,  extends  between  furmel  walls  and  is  parallel  to  the  X  axis.  The  gate  length,  {,  extends 
parallel  to  the  Y  axis;  thus,  groundwater  flows  parallel  to  the  Y  axis.  Finally,  Figure  5  suggests 
that  meters  were  used  as  the  length  dimension  in  the  FRAC-3D  simulations. 

Before,  the  numerical  model  could  be  formulated,  the  analytical  hydraulic  model  was 
used  to  configure  specific  funnel/gate  dimensions,  (  and  L,  using  the  following  input  data: 

Site  conditions: 

k  =  4.32  m/d 
(j),  =  5  m 

d(t)/dx  (X2  =  0)  =  -  0.005 

System  design  conditions: 

T=  ll.Od 
0  =  0.45 
k„  =  43.2  m/d 
B  =  61m 
b/B  =  0.6 
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Y  [m] 


Figure  5.  Conceptual  Orientation  of  the  Funnel/Gate  System  in  the  FRAC-3D 
Flow  Domain. 

Feasible  funnel/gate  dimensions  identified  with  the  analytical  model  were: 

()  =  0.9  m 
L  =  1.3  m 
b  =  36.6  m 

Keeping  Figure  5  in  mind,  the  groundwater  flow  simulation  from  FRAC-3D  can  be 
viewed  in  Figure  6.  This  figure  is  an  illustration  of  simulated  streamlines  created  in  the 
groundwater  flow  field  containing  the  above  funnel/gate  system  design.  The  numerical  values  on 
both  axes  are  distances  expressed  in  meters.  Groundwater  is  generally  flowing  parallel  with  the 
Y  axis.  The  funnel  walls  are  easily  visible  in  Figure  6,  and  the  gate  width  appears  equal  to 
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the  above  specified  design  width  of  36.6  meters.  The  FRAC-3D  simulation  results  show  the 
resulting  funnel/gate  design  intercepting  close  to  49  meters  of  flow  width  (as  viewed  from  the 
streamlines);  this  represents  80  percent  of  the  funnel  design  width  of  B  =  61  meters.  This 
suggests  that  the  anal)^ical  model  produced  a  funnel/gate  design  configuration  that  has  an  80 
percent  capture  efficiency.  It  also  says  the  volummetric  flow  rate  is  20  percent  less  than  the 
design  flow,  Q;  consequently,  the  numerically  simulated  hydraulic  residence  time  in  the  gate  will 
be  longer  than  the  critical  hydraulic  residence  time.  The  extended  residence  time  provides  for 
additional  contaminant  degradation;  thus,  the  system  design  is  conservative  in  the  regard  that  it 
will  produce  water  with  lower  contaminant  levels  than  the  design  quality.  These  results 
'demonstrate  that  the  analytical  hydraulic  model  can  identify  feasible,  albeit  conservative, 
fuimel/gate  designs  when  appropriate  information  is  provided. 

In  another  validation  effort,  a  second  funnel/gate  design  was  created  using  the  same  site 
and  system  conditions  as  described  above,  except  that  the  width  ratio  b/B  was  increased  to  0.7. 
The  result  of  increasing  the  gate  width,  b,  is  that  gate  length,  t,  will  decrease  such  that  the  total 
gate  volume  remains  unchanged;  thus,  the  total  volume  of  required  reactive  medium  does  not 
change.  This  happens  as  a  consequence  of  Equation  25  and  the  specified  critical  hydraulic 
residence  time  of  1 1 .0  days.  Under  the  new  width  ratio,  the  analytical  hydraulic  model  identified 
the  following  feasible  funnel/gate  dimensions: 

I  =  0.15  m 
L  =  .35m 
b  =  42.7  m 

Figure  7  is  an  illustration  of  the  groundwater  flow  streamlines  simulated  with  FRAC-3D 
for  a  second  funnel/gate  design.  Streamlines  show  this  funnel/gate  configiuration  intercepting  a 
groundwater  flow  width  of  52  meters  or  85  percent  of  funnel  design  width.  As  expected,  these 
results  show  that  increasing  the  gate  width  b  increases  the  capture  efficiency  of  the  system. 

Hence,  a  long  interception  trench  or  wide  permeable  wall  is  expected  to  be  more  efficient  than  a 
narrow  gate;  however,  this  says  nothing  about  comparative  costs  of  such  configurations.  These 
results  suggest  the  funnel/gate  design  model  could  be  used  to  pre-dimension  a  minimum  cost 
system;  however,  subsequent  numerical  simulations  should  be  conducted  to  verify  the  hydraulics 
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of  the  design.  If  these  numerical  simulations  demonstrate  that  the  proposed  design  is  unable  to 

intercept  a  plume  of  specified  width,  the  funnel  width.  B.  and/or  the  gate  width,  b,  is  increased  to 
enhance  capture. 

SUMMARY 

This  section  began  with  a  brief  review  of  an  existing  analytical  hydraulic  model  for 
funnel/gate  systems  developed  by  Christensen  and  Hatfield  (1 994).  From  this  an  iterative 
design  algonthm  was  developed  for  dimensioning  a  flinnel/gate  system.  The  design  algorithm 
incorporates  Equations  (24)  to  (28)  and  requires  the  following  site  and  system  data  to  execute: 

-  a)  the  critical  hydraulic  residence  time  in  the  gate; 

b)  the  transverse  vridth  of  groundwater  flow  to  be  intercepted  within  the  funnel,  B; 
this  is  practically  taken  to  be  greater  than  the  width  of  the  plume; 

c)  the  saturated  thickness  of  the  phreatic  aquifer,  (j),  at  the  entrance  of  the  funnel; 

d)  the  groundwater  flow  through  the  area  defined  by  width  B  and  flow  thickness,  (J), ; 

c)  the  aquifer  hydraulic  conductivity; 

f)  the  porosity  of  the  reactive  porous  media  within  the  gate;  and 

g)  the  hydraulic  conductivity  of  the  porous  media  within  the  gate. 

In  the  last  part  of  this  section,  results  of  two  FRAC-3D  numerical  validation  studies  were 
presented.  For  the  two  funnel/gate  systems  examined,  the  numerical  results  showed  the 
analytical  design  model  produced  system  configurations  which  exhibited  80  to  85  percent 
capture  efficiency.  These  efficiencies  increase  as  the  gate  becomes  wider.  Thus,  results  of  the 
numencal  validation  suggest  the  funnel/gate  design  model  could  be  used  to  pre-dimension  a 

system;  however,  subsequent  numerical  simulations  should  be  conducted  to  verify  the  hydraulics 
of  the  design. 
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SECTION  IV 


COST  MINIMIZATION  MODELING 

The  primary  objective  of  this  chapter  is  to  develop  a  cost  estimation  model  that  will 
identify  minimum  cost  funnel/gate  designs.  The  secondary  objective  is  to  recast  the  general  cost 
model  into  Lagrangian  optimization  formulation  that  is  then  amenable  to  solution. 

COST  OBJECTIVE  FUNCTION 

The  total  cost  of  funnel/gate  systems  can  include  at  least  the  five  components: 

1)  Construction/design  costs  per  unit  length  of  funnel  wing  walls; 

2)  Construction/design  costs  per  unit  length  of  gate; 

3)  Construction/design  costs  per  unit  Width  of  gate; 

4)  Land  costs  per  unit  area;  and 

5)  Costs  of  reactive  media  installed  in  the  gate. 

Additional  costs  such  as  water  quality  monitoring  before  and  after  system  installation  are  not 
presently  considered. 

Construction/design  cost  per  imit  length  of  funnel  wing  walls  can  be  summarized  as 
follows: 

Total  cost  of  4  funnel  walls  =  (29) 


where  Xj  =  the  length  of  a  single  fuimel  wall;  and  E,  =  the  cost  per  unit  length  of  funnel  wall. 
This  unit  cost  is  calculated  as  follows: 

d  (30) 


in  which  d  =  depth  of  funnel  walls  (i.e.  the  depth  to  which  slurry  walls  are  constructed  or  the 
depth  to  which  sheet  piles  are  driven);  and  C^ffw  =  cost  per  unit  area  of  sheet  pile  or  slurry  wall 


24 


used  to  construct  the  funnel  walls.  Clearly  the  cost  parameter,  ,  varies  with  the  total  area  of 
sheet  pile  used  and  the  depth  to  which  the  sheet  pile  is  driven.  Table  1  presents  Csff«,  values  for 
sheet  pile  where  costs  vary  as  a  function  of  both  total  sheet  pile  area  installed  and  pile  depth.^ 


TABLE  1.  VALUES  FOR  SHEET  PILE. 


Sheet  Pile  Unit  Costs  Csfi^  [  Dollars  ] 

Sheet  Pile 

Sheet  Pile  Area  [  m^  ] 

Depth  [  m  ] 

<465 

465  -  860 

>  1860 

0-8 

280 

193 

172 

8-15 

301 

215 

183 

15-24 

322 

236 

204 

Construction/design  costs  per  unit  length  of  gate  are  quantified  as  follows; 

Total  cost  of  2  longitudinal  gate  walls  =  2  (31) 


where  H  =  the  longitudinal  gate  length  as  defined  in  Section  3,  and  Ej  =  the  cost  per  unit  length 
of  longitudinal  gate  wall.  E2  estimated  from: 

P2) 


in  which  Cjfgi  =  cost  per  unit  area  of  sheet  pile  or  slurry  wall  used  to  construct  the  longitudinal 


^Jowett,  R.,  Personal  Communication,  Waterloo  Educational  Services,  Inc.,  2  Taggart 
Court,  Unit  4,  Guelph,  Ontario  NIH  6H8, 1996. 
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gate  walls.  Thus,  £3  may  equal  E,  if  gate  and  funnel  walls  are  constructed  using  similar 
materials,  i.e.,  =  C^fg,. 

Construetion/design  cost  per  unit  width  of  the  transverse  gate  wall  can  be  summarized  as 
follows: 

Total  cost  of  2  transverse  gate  walls  =  2E^BX^  (33) 


where  X,  =  b/B,  the  dimensionless  ratio  of  gate  width  to  funnel  width;  B  =  the  funnel  wall 
entrance  width  as  defined  in  Section  3;  and  £3  =  the  cost  per  unit  width  of  the  transverse  gate 
wall.  This  cost  parameter  is  estimated  from: 


^3 


(34) 


in  which  C^fg^  =  cost  per  unit  area  of  sheet  pile  or  trenching  required  to  construct  the  transverse 
gate  wall. 

Land  costs  are  estimated  with: 

Total  Land  Costs  =  +  2L)  (35) 

where  L  =  the  projected  length  of  the  funnel  walls  as  defined  in  Section  3;  and  En  =  land  costs 
per  unit  area. 

Finally,  there  is  the  cost  of  the  reactive  medium  installed  in  the  gate.  This  cost  is 
essentially  constant  within  a  suite  of  alternative  designs  being  examined  for  any  single  site.  This 
is  because  the  critical  hydraulic  residence  time  must  be  preserved  between  the  different  designs 
under  investigation;  it  may  be  recalled  from  Section  2  that  to  allow  the  hydraulic  residence  time 
to  vary  between  potential  designs  would  produce  different  systems  which  are  not  comparable 
because  each  will  achieve  a  different  degree  of  contaminant  degradation.  In  keeping  with  a 
constant  critical  hydraulic  residence  time,  the  gate  volume  becomes  a  known  constant  (so  long 
as  Q  remains  constant).  The  cost  of  the  reactive  gate  material  is  estimated  from: 
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^12  07’ 

Total  cost  of  gate  reactive  media  =  — - — 


(36) 


in  which  Q  =  the  volumetric  flow  through  the  gate  as  defined  in  Section  3,  T  =  critical  hydraulic 
residence  time  defined  in  Sections  2  and  3,  and  £,2  =  cost  of  reactive  media  per  unit  volume. 

The  cost  of  iron  is  generally  $1,690  -  1,760  m'^  (Jowett).^ 

Equations  (29),  (31),  (33),  (35),  and  (36)  ean  be  combined  to  obtain  an  overall  cost 

objective  function. 

Total  funnellgate  costs  =  4E^X2  +  lE-^BX^  +  E^^B{ll  +  2L)  +  E^2QT  (37) 

CONSTRAINTS 

The  minimum  value  of  the  cost  objective  function.  Equation  (37),  is  determined  imder  a 
set  of  constraints  which  ensure  that  only  feasible  funnel/gate  designs  are  identified  (a  feasible 
design  being  one  that  is  correct  both  from  hydraulic  and  from  transport  perspectives).  In  this 
study,  four  constraints  were  considered: 

1)  Constraint  identifying  feasible  combinations  of  gate  dimensions  length,  {,  width  b, 
and  saturated  thickness,  (Jjgatej  such  that  the  critical  hydraulic  residence  time  is  achieved  inside  the 
gate; 

2)  Constraint  to  identify  feasible  gate  potentials  (saturated  thickness),  (Jjgates  that  are 
eonsistent  with  the  hydraulic  model; 

3)  Constraint  defining  funnel  walls  lengths  in  terms  of  gate  width,  b,  and  projected 
funnel  wall  length,  L; 

4)  Constraint  to  identify  feasible  combinations  of  projected  funnel  wall  length,  L, 
gate  width,  b,  and  gate  length, «,  consistent  with  the  hydraulic  model. 


^Jowett,  R.,  Personal  Commumcation,  Waterloo  Educational  Services,  Inc.,  2  Taggart 
Court,  Unit  4,  Guelph,  Ontario,  NIH  6H8, 1996. 
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The  first  constraint  functions  to  ensure  feasible  combinations  of  gate  dimensions  width, 
b,  length, and  saturated  thickness,  4)ga,e.  are  identified  to  produce  a  saturated  gate  volume  large 
enough  to  meet  the  critical  hydraulic  residence  time,  T.  It  may  be  recalled  that  the  value  of  T  is 
obtained  from  executing  the  transport  model.  It  represents  the  longest  "adequate  contact  time"m 
the  gate  which  ensures  all  target  contaminants  meet  desired  water  goals  at  the  gate  exit. 

Recalling  Equation  (25): 


^  gate 


Q 


=  T 


(38) 


or 


5  0 


(39) 


where 


-  ^gate  (40) 

Constraint  Equation  (39)  ensures  that  the  gate  volume  will  be  large  enough  to  produce  a 
hydraulic  residence  time  equal  to  the  specified  critical  hydraulic  residence  time. 

The  second  constraint  is  required  to  ensure  the  cost  model  identifies  feasible  potentials 
(saturated  thickness)  in  the  gate  that  are  consistent  with  the  hydraulic  model.  This  constraint  is 
derived  from  Equation  (28): 


2QL 

k(B  -  b) 


(41) 


Rearranging  Equation  (41)  gives  Equation  (42),  the  final  form  of  the  second  constraint  on  the 
gate  potential: 
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^/(i  -  X,)  -  (1),2(1  -  X,)  -  =  0. 


(42) 


The  cost  objective  function,  Equation  (37),  is  defined  in  terms  of  variable  X2,  the  fimnel 
wall  length;  consequently,  the  third  constraint  is  required  because  it  defines  Xj  in  terms  of  design 
variables,  projected  funnel  wall  length,  L,  and  the  dimesionless  width  ratio,  X,.  Thus,  fi-om  the 
funnel/gate  system  geometry  the  following  constraint  was  derived: 

-  ^(1  -  2X^  +  =  0.  (43) 


The  fourth  and  final  constraint  is  derived  from  Equation  (26)  rewritten  below  to  include 
decision  variable,  Xi  : 


^  1 

L  -  1-^1— 
2  \n(X,)  ^  ^ 

1  -  (^1)  ^ 


(44) 


Equation  (44)  serves  to  limit  the  optimization  model  to  examine  only  those  combinations  of 
projected  funnel  wall  length,  L,  gate  width,  b,  and  gate  length,  f,  that  will  also  meet  hydraulic 
constraints  defined  by  governing  flow  equations  (i.e.,  continuity  and  momentum).  A  more 
convenient  form  for  Equation  (44)  is: 

2X,iX,ln[X^]  +  X,  -  X,2)  -  X3(A[1  -  -  X,  +  X,^)  =  0.  (45) 


Minimizing  the  objective  function  Equation  (37),  subject  to  the  fom  constraint  equations. 
Equations  (39),  (42),  (43),  and  (45),  constitutes  the  funnel/gate  minimum  cost  design 
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optimization  problem.  The  complete  nonlinear  formulation  is  summarized  below  and  contains 
five  primary  decision  variables,  Xj,  X2,  S,  L,  and  X5. 

Total  funnel! gate  costs  =  4E^X2  +2^2®  +  £jj5(C  +  2L)  +  E^2QT  (46) 


subject  to: 

BXfX^  -  ^  =  0 
*  ^  0 


(47) 


-X.)-$,2(l  -  X,)  -  ^\n[X,]  -  0. 


(48) 


X/  -  -  2X,  +  X,^)  -  L- 


0. 


(49) 


2X,{Xfn{X,-\  +  X,  -  -  JT3 


f[l  -X,]  -X,  -X,^ 


=  0. 


(50) 


LAGRANGIAN  OPTIMIZATION 

The  above  cost  minimization  formulation  is  a  nonlinear  optimization  problem.  Because 
this  formulation  contains  only  equality  constraints  (i.e.,  no  inequalities),  a  solution  can  be 
obtained  readily  by  recasting  the  problem  as  a  Lagrangian  optimization  formulation.  The  first 
step  of  the  approach  is  to  formulate  the  Lagrangian,  L^,  by  combining  the  objective  function  and 
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all  the  constraints.  Thus, 


AE^X^  +2E^ll  +  lE^BX^  +  E^^B(ll  +  21)  +  E^^QT 


+  X 


BX,IX,  -  ^ 


tX, 


X,\l  -X,)  -<1>,2(1  -X,)  - 


2QL 

kB 


In[X] 


+  X 


\X^2  -  i-{\  -  2X,  +X,2)  -  I- 


+  X 


2X5(Xiln[XJ  +  JTj 


-X3(A[1  -X,]  -X,  +X,2) 


(51) 


where  X^,  Xj,  Xg,  and  X9  =  Lagrangian  multipliers. 

The  next  step  to  reformulating  the  optimization  problem  is  to  derive  nine  equations  from 
differentiating  L,  with  respect  to  each  of  the  five  original  decision  variables  and  each  of  the  four 
Lagrangian  multipliers. 


dL 


F,  =  ^  =  2E.B  +  X^BHX^  -  X^ 


dX, 


V 2  _  ^  2  +  2QL 
5  1  uipv 


kBX^ 


X^B‘ 


(2  -  IX,) 
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2LX,(\b[X,]  +  2.  -  2X,)  +  IX,  ^  -  2X,  =  0 


WJ-/ 

^2  =  Trf  =  ^.  +  2XX  =  0 


'  if  =  *  ^nB  *  BX,X,X,  -  X,  -fd  -  X,)  -  X,  +  X,^  =  0 


^4  =  ^  =  2£„S  -  2iXj  4  2X,[X,ln(X,)  4  X,  -  X^j  =  0 


^6  =  ^  =  SX,X,I  -  Ig  = 


^7  "  —  = 


-  Oj^xi  -  Z,)  -  ^^ln{X^)  =  0 

kB 
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(59) 


^  T  2 

^  =  X"  -  —(1  -  2X  +  X^)  -£2=0 
*  aXg  ^  4  '  ‘ 


dL„ 

X  =  — ^  =  2Z 

^  ax 


X,ln(X,)  +  -  X,' 


f(l  -Z,)  -X, 

K 


=  0 


(60) 


Equations  (52)  -  (60)  constitute  the  Lagrangian  cost  optimization  model  comprised  of 
nine  nonlinear  equations  with  nine  unknowns  (the  five  original  decision  variables  and  the  four 
Lagrangian  multipliers).  These  equations  can  be  solved  using  a  generalized  Newton-Raphson 
algorithm  (Smith,  Henton,  and  Lewis,  1983).  Assuming  a  solution  is  found  to  the  above  nine 
equations,  that  solution  represents  a  stationary  point  that  may  or  may  not  represent  the  minimum 
cost  solution.  To  facilitate  the  Newton-Raphson  algorithm  in  the  search  for  the  minimum  cost 
solution,  initial  estimates  of  all  variables  are  chosen  which  are  thought  to  reflect  values  close  to 
the  minimum  cost  solution.  Hence,  the  Newton-Raphson  algorithm  begins  an  iterative  search  for 
the  optimum  funnel/gate  design  (lowest  cost)  beginning  with  variable  values  that  reflect  a  close 
approximation  of  optimum  solution. 

SUMMARY 

In  the  first  half  of  this  Section  a  funnel/gate  cost  estimation  model  was  developed  to 
identify  minimum  cost  designs  subject  to  four  constraints:  one  that  defines  the  desired  hydraulic 
retention  time  in  the  gate,  one  that  defines  the  funnel  wall  length  in  terms  of  other  system 
dimensions,  and  two  that  serve  as  hydraulic  constraints.  Costs  considered  in  the  model  include 
land  costs,  costs  per  unit  length  of  funnel  walls,  costs  per  unit  length  of  gate,  costs  per  unit  width 
of  gate,  and  costs  associated  with  the  reactive  medium  used  in  the  gate. 

The  second  half  of  this  Section  was  devoted  to  recasting  the  general  cost  model  into  a 
Lagrangian  optimization  model.  As  a  result  of  this  reformulation,  the  optimization  problem  was 
reduced  to  that  of  solving  a  system  of  nine  nonlinear  equations  [Equations  (52)  to  (60)]  with 
nine  unknowns.  Recommendations  were  given  as  to  how  the  nine  equations  could  be  solved. 
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SECTION  V 


funnel/gate  design  model  (FGDM) 

In  this  chapter  two  objectives  will  be  achieved.  The  first,  is  to  provide  a  brief  overview 

of  the  Funnel/Gate  Design  Model  (FGDM)  software.  The  second  objective  is  to  step  through  an 
example  problem  using  FGDM. 

FGDM 


FGDM  is  a  FORTRAN  program  developed  to  ascertain  feasible  dimensions  of 
funnel/gate  systems  and  also  identify  the  lowest  cost  design.  Feasible  is  defined  here  to  include 
designs  that  are  consistent  with  the  transport  and  hydraulic  theories  outlined  in  Sections  2  and  3. 
The  software  includes  a  primary  control  program,  FGDM.for,  and  four  subroutines,  TRANS.for, 
HYD.for,  OPTI.for  and  GAUSSJ.for.  During  execution,  subroutine  TRANS.for  is  first  called, 
then  HYD.for,  and  finally  OPTI.for.  GAUSSJ.for  is  a  Gauss-Jordan  elimination  algorithm  with 
full  pivoting  that  was  developed  by  Press,  Flannery,  Teukolsky,  and  Vetterling  (1990): 
GAUSSJ.for  is  called  directly  by  OPTI.for.  Each  subroutine  writes  results  to  a  single  output 
file.  The  structure  of  the  program  is  such  that  new  subroutines  may  be  substituted  for  existing 
versions  without  affecting  other  components  of  FGDM.  The  program  requires  one  input  file  and 
produces  one  output  file.  In  addition,  it  writes  to  the  screen  much  of  the  same  design  information 
written  to  the  output  file.  Section  6  provides  a  complete  listing  of  the  program  along  with 
examples  of  input  and  output  files. 

FGDM.for  is  the  primary  control  program  that  reads  the  input  file,  dimensions  several 
common  block  arrays  and  calls  subroutines  TRANS.for,  HYD.for,  and  OPTI.for.  This  program 
also  prompts  the  user  for  both  input  and  output  file  names.  Initially  FGDM.for  opens  and  reads 
the  input  file  and  then  calls  subroutine  TRANS.for  to  perform  contaminant  fate  and  transport 
simulations  in  the  gate.  TRANS.for  returns  control  to  FGMD.for  with  an  estimate  of  the  critical 
hydraulic  residence  time  for  the  gate.  Next,  FGMD.for  calls  the  hydraulic  design  subroutine 
HYD.for  to  develop  several  funnel/gate  designs  using  the  aforementioned  hydraulic  residence 
time  in  each  design.  From  a  list  of  width  ratios,  b/B,  specified  in  the  input  file,  HYD.for 
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formulates  one  system  design  for  each  width  ratio.  HYD.for  also  estimates  the  cost  of  each 
design  and  then  returns  control  to  FGMD.for.  Finally,  FGMD.for  calls  subroutine  OPTI.for  to 
search  for  the  minimum  cost  system.  OPTI.for  initiates  this  search  using  the  dimensions  of  the 
lowest  cost  system  identified  by  HYD.for.  Once  OPTI.for  finds  the  true  minimum  cost  system 
configuration,  it  returns  control  to  FGDM.for  and  the  program  terminates.  Additional  details 
regarding  FGDM.for  can  be  obtained  reading  the  documented  code  (See  Section  6). 

TRANS.for  is  a  contaminant  fate  and  transport  subroutine  that  ascertains  the  critical 
hydraulic  residence  time,  T,  in  the  gate  using  an  iterative  algorithm.  The  subroutine  begins  the 
iterative  search  by  equating  the  hydraulic  residence  time,  T  =  DT,  where  DT  =  a  time  step 
specified  in  the  input  file.  Next,  the  TRANS.for  solves  Equations  (5),  (8),  (1 1),  and  (17)  and 
then  checks  to  determine  whether  calculated  values  for  C„  Cj,  C3,  or  C4  exceed  respective  water 
quality  standards  STDl,  STD2,  STDS,  or  STD4.  If  any  of  the  four  standards  are  exceeded,  the 
subroutine  increases  the  value  of  T  by  an  increment  equal  to  DT  (now  T  =  2-DT)  and  then 
reevaluates  Equations  (5),  (8),  (1 1),  and  (17)  for  comparison  against  water  quality  standards. 

This  iterative  process  continues  (and  the  value  of  T  is  increased  after  each  iteration)  so  long  as 
any  of  the  calculated  contaminant  levels,  C„  Cj,  C3,  or  C4  exceed  their  respective  water  quality 
standard.  The  iterative  search  stops  when  a  value  T  is  found  that  produces  Cj,  Cj,  C3,  and  C4 
values  less-than-or-equal-to  pertinent  water  quality  standards.  When  the  search  has  terminated, 

T  =  the  critical  hydraulic  residence  time.  During  the  search  for  this  residence  time,  TRANS.for 
writes  to  an  output  file  C,,  C2,  C3,  and  C4  values  calculated  for  each  value  of  T  evaluated. 

Finally,  before  it  returns  control  to  FGFM.for,  it  writes  to  the  output  file  the  value  of  the  critical 
hydraulic  residence  time. 

HYD.for  is  the  hydraulic  subroutine  that  generates  feasible  funnel/gate  designs.  In  this 
application,  the  subroutine  takes  the  critical  hydraulic  residence  time  produced  by  TRANS.for  to 
generate  several  systems  designs.  Specifically,  the  subroutine  solves  Equations  (24)  through  (28) 
using  the  iterative  design  algorithm  described  in  Section  3.  Each  design  is  created  using  a  unique 
width  ratio,  b/B,  obtained  fi’om  a  list  specified  in  the  input  file.  From  the  list  of  specified  b/B 
values,  only  ratios  greater  than  the  hydraulic  conductivity  ratio,  k/k,,,  are  used;  it  may  be  recalled 
from  Christensen  and  Hatfield  (1994)  that  infeasible  designs  are  obtained  for  b/B  <  k/kf,. 
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HYD.for  also  calculates  the  cost  of  each  design  using  Equation  (37).  During  execution,  HYD.for 
writes  to  the  output  file  the  cost  and  design  dimensions  of  each  fiinnel/gate  system  configured. 

OPTI.for  is  the  subroutine  used  to  identify  the  minimum  cost  funnel/gate  design. 
Essentially  this  subroutine  employs  a  generalized  Newton-Raphson  algorithm  (Smith,  Hinton, 
and  Lewis,  1983)  combined  with  a  Gauss-Jordan  elimination  algorithm  with  full  pivoting  (Press, 
Flannery,  Teukolsky,  and  Vetterling,  1990)  to  solve  the  system  of  nine  nonlinear  equations 
[Equations  (52)  -  (60)  ]  of  the  Lagrangian  cost  model.  OPTI.for  initiates  an  iterative  search  for 
the  minimum  cost  design  using  the  dimensions  of  the  lowest  cost  system  identified  by  HYD.for. 
During  the  search  process,  OPTI.for  calls  subroutine  GAUSS  J.for  to  solve  a  new  system  of 
equations  generated  from  the  application  of  the  Newton-Raphson  algorithm  to  the  original 
aforementioned  nine  nonlinear  equations.  Once  OPTI.for  finds  the  true  minimum  cost  system 
configuration,  the  cost  and  dimensions  of  this  system  are  written  to  the  screen  and  to  the  output 
file.  Finally,  the  subroutine  writes  to  the  output  file  the  surface  areas  associated  with  the  funnel 
walls,  the  longitudinal  gate  walls,  and  the  transverse  gate  walls.  These  surface  areas  can  be  used 
to  verify  that  appropriate  values  have  been  used  to  define  cost  parameters,  Cjfgi,  Cshv^,  and  Csfg^,  in 
the  input  file. 

EXAMPLE  PROBLEM 

To  demonstrate  FGDM,  a  problem  example  was  formulated.  Consider  a  site  with  a 
plume  that  is  49  meters  wide.  The  plume  contains  PCE  and  TCE  with  a  maximum  concentration 
of  each  at  5000  pg/1.  The  local  phreatic  surface  is  0.5  meters  below  ground  surface,  and  the 
depth  to  the  aquiclude  is  5.5  meters;  thus,  (j)]  =  5  m  .  The  local  groundwater  flow  gradient  is 
d(j)/dx  =  -0.005  in  an  aquifer  with  an  estimated  hydraulic  conductivity  of  k  =  4.32  m/d. 

The  suggested  groundwater  remediation  approach  for  the  site  is  to  intercept  and  treat  the 
polluted  groundwater  using  a  funnel/gate.  The  applicable  system  performance  standard  is  for 
water  exiting  the  gate  to  meet  water  quality  standards  of  5  pg/1  PCE,  5  pg/1  TCE,  70  pg/l  DCE, 
and  1  pg/1  vinyl  chloride. 

To  initiate  an  investigation  of  feasible  system  designs,  typical  funnel/gate  parameters  are 
assumed  (i.e.,  degradation  parameters,  gate  porosities,  gate  widths,  and  gate  hydraulic 
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conductivities).  The  gate  porosity  and  hydraulic  conductivity  are  assumed  to  be,  respectively,  6 
=  0.45  and  k,,  =  43.2  m/d.  The  gate  width,  b,  is  determined  by  the  value  of  the  funnel  width  B 
and  the  gate  to  funnel  width  ratio,  b/B.  FGDM  can  be  used  to  study  a  suite  of  up  to  10 
funnel/gate  designs.  Each  design  reflects  a  unique  b/B  ratio;  however,  FGDM  ignores  ratios  with 
values  less  than  the  hydraulic  conductivity  ratio,  k/ko.  As  for  the  funnel  width,  B,  this  parameter 
is  estimated  in  advance  assuming  a  reasonable  capture  efficiency,  i.e.,  80  percent  (see  example 
validation  problems  presented  in  Section  3).  Under  this  assumption,  a  funnel  entrance  width  B  = 
61  meters  is  needed  to  intercept  a  plume  49  meters  wide.  Finally,  there  are  degradation 
parameters,  Gillham  (1995)  suggests  the  following  first-order  estimates:  =  36.9  d"'  and  A,,b  = 
0.0  d-'  for  PCE;  X2,  =  3.69  d*'  and  Xjb  =  33.18  d’’  for  TCE;  A3,  =  .83  d*'  and  X^^  =  0.0  d''  for 
PCE;  and  A,,  =  .83  d"'  for  vinyl  chloride. 

In  this  example,  it  is  assumed  that  both  funnel  and  gate  structures  will  be  constructed 
from  sheet  piles.  Since  the  aquiclude  exists  5.5  meters  below  ground  surface,  these  piles  will  be 
driven  to  a  total  depth,  d  =  6  meters  ( allowing  0.5  meters  of  sheet  pile  depth  as  a  key  into  the 
aquitard).  The  values  of  applicable  sheet  pile  cost  parameters  can  be  obtained  from  Section  4. 
The  reactive  gate  material  costs  are  taken  to  equal  $1,725  m  while  land  costs  are  assumed 
small. 

The  above  information  was  incorporated  into  an  FGDM  input  file  shown  on  the  next 
page.  Line  numbers  were  added  to  the  right  of  the  file  text  for  quick  referencing.  Line  3 
contains  water  quality  standards  STDl  (5  |ag/l),  STD2  (5  pg/1),  STD3  (70  |ag/l),  and  STD4  (1 
pg/1)  for  respective  contaminants  Cj  (PCE),  Cj  (TCE),  C3  (DCE),  and  C4  (vinyl  chloride).  The 
last  number  appearing  in  line  3  is  the  time  step,  DT  (0. 1  days),  used  in  the  subroutine 
TRANS.for,  to  increase  the  hydraulic  residence  time  between  iterative  calculations  of  C,,  C2,  C3, 
and  C4.  The  assumed  funnel  width,  B  (61  m),  the  hydraulic  gradient,  d4)/dx  (-0.005),  the 
saturated  aquifer  thickness,  (J)i  (5  m),  and  the  aquifer  hydraulic  conductivity,  k  (4.32  m/d)  are 
written  in  line  7.  The  gate  porosity,  6  (0.45),  the  aquifer  gate  hydraulic  conductivity  ratio  k/k,, 
(0.1),  and  the  depth  of  the  funnel/gate  walls  d  (6  meters)  are  listed  on  line  11.  Lines  17  and  21 
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Example  Input  File: 

Data 

STDl  [ug/1]  STD2  [ug/1]  STDS  [ug/1]  STD4  [ug/1]  DT(TIMESTEP)  PAYS] 

* 

5.0  5.0  70.0  1.0  .1 

* 

FUNNEL  WIDTH  [M]  GRADIENT  AQU.  THICKNESS  [M]  AQU.  HYDRAULIC  COND.  [M/d] 

* 

61.  -0.005  5.  4.32 

* 

GATE  POROSITY  HYDRAULIC  COND.  RATIO  (AQUIFER/GATE)  DEPTH  OF  SYSTEM  WALLS  [M] 

* 

.45  .1  6. 

* 

VARIOUS  RATIOS  OF  GATE  WIDTH  TO  PLUME  WIDTH  b/B  EXAMINED 
* 

BR(1)  BR(2)  BR{3)  BR(4)  BR(5) 

* 

.02  .05  .1  .2  .4 


Line  Number 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


BR(6)  BR(7)  BR(8)  BR(9)  BR(IO) 

♦ 

.5  .6  .7  .8  .9 

♦ 

♦♦♦♦DATA  PLUME  CONTAMINANT  CONCENTRATIONS  AND  DEGRADATION  PARAMETERS 


PLUME  Cl  CONC.  [ug/1] 

* 

LAMBDA  lA  [1/DAY] 

LAMBDA  IB  [1/DAY] 

5000. 

* 

36.9 

0. 

PLUME  C2  CONC.  [ug/1] 

LAMBDA  2A  [1/DAY] 

LAMBDA  2B[1 /DAY] 

5000. 

* 

3.69 

33,18 

PLUME  C3  CONC.  [ug/1] 

♦ 

LAMBDA  3A  [1/DAY] 

LAMBDA  3B  [1/DAY] 

0. 

.8318 

0. 

PLUME  C4  CONC.  [ug/1] 

* 

LAMBDA  4A[I/DAY] 

0. 

* 

.8318 

♦  READ  IN  COSTS  COEFFICIENTS:  [  IN  UNITS  OF  THOUSANDS  OF  DOLLARS  ] 

♦  CSFFW  =  COST  PER  UNIT  AREA  OF  FUNNEL  WALL  MATERIAL 

♦  CSFGL  -  COST  PER  UNIT  AREA  OF  LONGITUDINAL  GATE  WALL 

♦  CSFGW  =  COST  PER  UNIT  AREA  OF  TRANSVERSE  GATE  WALL 

♦  E 1 1  =  LAND  COSTS  PER  UNIT  AREA 

♦  E 1 2  =  REACTIVE  MEDIA  COSTS  PER  UNIT  VOLUME 


19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 


♦  47 

CSFFW  CSFGL  CSFGW  Ell  E12  48 

*  49 

.193  .193  .193  .001  1.725  50 

♦  51 

*  READ  IN  INITIAL  ESTIMATES  OF  LAGRANGIAN  MULTIPLIERS,  X6,  X7,  X8,  AND  X9  52 

*  53 

X6  XI  X8  X9  54 

♦  55 

1.  1.  1.  10.  56 
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comprise  a  listing  of  width  ratios,  b/B,  to  be  evaluated  with  FGDM.  The  maximum  PCE 
concentration  C,o  =  5000  pg/1  and  pertinent  decay  parameters  =  36.9  d'*  and  X,b  =  0.0  d''  for 
PCE  are  written  on  line  27.  On  line  3 1  there  is  the  maximum  TCE  concentration  C20  =  5000 
pg/l  and  associated  decay  parameters  ^23  =  3.69  d"'  and  X2\>  =  33.18  d '.  Line  35  lists  the 
maximum  plume  concentration  for  DCE  C30  =  0  pg/l  3nd  decay  parameters  —  .83  d  and  A,3(j 
=  0.  d  '.  For  vinyl  chloride,  input  file  line  39  contains  the  maximum  plume  concentration  C40  = 

0  pg/1  and  decay  parameter  X^,  =  .83  d  ' .  The  next  important  data  line  is  number  50.  This  line 
contains  values  for  cost  parameters  Csflvv  =  -193,  C^fg,  =  .193,  and  C^fg^  =  .193.  In  addition,  this 
line  includes  the  cost  per  unit  area  of  land,  E„  =  .001,  and  the  cost  per  unit  volume  of  reactive 
gate  material,  E,2  =  1.725  thousand  dollars/m^  Note,  all  cost  coefficients  are  expressed  in  terms 
of  thousands  of  dollars  (compare  with  Table  1).  Finally,  there  is  line  56  which  contains  initial 
estimates  of  the  values  of  Lagrangian  multiplier  =1 .,  X7  =  1 .,  Xg  =  1 .,  and  X9  =  10.  Unless 
there  are  severe  convergence  problems,  line  56  need  not  be  changed  between  applications.  If 
these  problems  do  occur,  it  is  suggested  that  the  user  look  at  the  output  file  first.  Output  will 
exist  from  HYD.for,  which  may  reveal  that  a  low  cost  system  doesn't  exist.  Making  a  plot  of 
system  eosts  versus  system  width  ratios,  b/B,  (see  Figure  10  next  section)  may  reveal  a  flat  or 
monotonic  curve  that  has  no  minimum  value;  this  would  explain  why  OPTI.for  failed  to 
converge  on  a  low  cost  system. 

Using  the  above  input  file,  FGDM  generates  an  output  file  that  may  be  divided  into  the 
three  respective  sections  produced  from  subroutines  TRANS.for,  HYD.for  and  OPTI.for.  Shown 
on  the  next  page  is  an  abbreviated  output  file  from  the  above  example  problem.  The  first  section 
is  generated  from  subroutine  TRANS.for.  This  section  contains  five  printed  columns,  the  first 
being  the  hydraulic  residence  time  T,  and  the  remaining  four  reflecting  calculated  values  of 
contaminants  C„  C2,  C3,  and  C4  ( in  this  case  PCE,  TCE,  DCE,  and  vinyl  chloride).  The  end  of 
the  first  section  is  marked  by  a  statement: 

CRITICAL  HYDRAULIC  RESIDENCE  TIME  =  11.10  [DAY] 

This  statement  indicates  the  value  of  the  critical  hydraulic  residence  time  in  the  gate  determined 
by  the  conditions  defined  in  the  example  problem.  Section  one  results  showed  PCE,  TCE,  and 
DCE  concentrations  decreased  to  performance  standards  in  the  gate  long  before  vinyl  chloride; 
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however,  by  day  11.1  vinyl  chloride  decreased  below  the  1 .0  |ig/l  standard.  Figxne  8  is  a  plot  of 
the  data  from  the  first  section  showing  PCE,  TCE,  DCE,  and  vinyl  chloride  concentrations  as  a 
flmction  of  gate  hydraulic  residence  time.  Figure  8  clearly  shows  the  degradation  of  PCE  and 
TCE  and  the  subsequent  production/degradation  of  DCE  and  vinyl  chloride.  Figure  9  contains 
only  the  first  two  days  of  data  used  in  Figure  8.  This  figure  illustrates  more  clearly  the  rapid 
disappearance  of  PCE  and  TCE. 

The  second  section  of  the  output  file  is  generated  by  subroutine  HYD.for.  Each  line  of 
this  section  contains  the  dimensions  and  the  cost  of  a  single  feasible  funnel/gate  system  for  a  b/B 
ratio  specified  in  the  input  file.  It  should  be  noted  that  b/B  ratios  0.02,  0.05,  and  0.1  were  listed 
in  the  input  file,  but  output  lines  were  not  generated  by  HYD.for  because  they  are  less  than 
k/ko  =  0.1 .  The  first  three  columns  of  output  in  this  section  are  gate  widths,  b,  lengths,  i,  and  gate 
volumes.  The  basic  length  unit  is  meters.  It  may  be  observed  that  the  gate  volume  is  constant 
between  designs,  while  gate  length  and  width  vary.  A  constant  gate  volume  is  expected  because 
the  hydraulic  residence  time  is  a  constant.  The  next  two  columns  give  fiinnel  wall  lengths,  Xj 
and  projected  funnel  wall  lengths,  L.  Under  the  next  heading  of  ‘Dimensionless  Parameters’ 
there  is  a  column  for  gate/fiuinel  width  ratios,  b/B  and  another  for  length  ratios,  L/C;  these 
parameters  found  relevance  in  Christensen  and  Hatfield  (1994).  Finally,  the  last  column  contains 
total  system  costs  expressed  as  thousands  of  dollars.  Figure  10  is  a  plot  of  the  total  system  cost 
versus  the  width  ratio,  b/B.  Clearly,  Figure  10  identifies  the  lower  cost  systems  as  those  with 
gate/funnel  width  ratios  ranging  from  0.5  to  0.7.  Given  that  funnel/gate  capture  efficiency 
increases  with  the  value  of  the  width  ratio,  b/B  (recall  validation  studies  in  Section  3),  it  may  be 
desirable  to  select  among  the  lower  cost  designs,  the  one  with  the  largest  width  ratio. 

The  beginning  of  the  third  and  last  section  of  the  output  file  is  marked  by  the  heading: 

*********  MINIMUM  COST  DESIGN  ******** 

This  section  is  written  by  subroutine  OPTI.for  and  has  initially  the  same  format  as  the  second 
output  section.  The  third  section  initially  lists  the  dimensions  and  cost  of  the  true  minimum  cost 
design.  The  total  cost  of  this  design  is  $424,000,  and  it  has  the  following  dimensions: 
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Example  Output  File: 


TIME 

Cl 

C2 

C3 

C4 

[DAY] 

[ug/I] 

[ug/1] 

[ug/1] 

[ug/1] 

.00 

5000. 

5000. 

.00 

.00 

.10 

124.8 

586.6 

879.91 

48.15 

.20 

3.12 

26.22 

873.98 

115.50 

.30 

.08 

.94 

806.89 

173.34 

.40 

.00 

.03 

742.58 

221.27 

.50 

.00 

.00 

683.31 

260.45 

.60 

.00 

.00 

628.77 

291.96 

.70 

.00 

.00 

578.59 

316.79 

.80 

.00 

.00 

532.41 

335.79 

10.70 

.00 

.00 

.14 

1.25 

10.80 

.00 

.00 

.13 

1.16 

10.90 

.00 

.00 

.12 

1.08 

11.00 

.00 

.00 

.11 

1.00 

11.10 

.00 

.00 

.10 

.93 

CRITICAL  HYDRAULIC  RESIDENCE  TIME  =  1 1 . 1 0  [  DAY  ] 


GATE 

FUNNEL 

DIMENSIONLESS 

TOTAL  COSTS 

PARAMETERS 

WIDTH 

LENGTH 

VOLUME 

LENGTH 

PROJECTED 

b/B 

L/1 

LENGTH 

THOUSANDS  OF 

[M] 

[M] 

[M^3] 

[M] 

[M] 

DOLLARS 

12.20 

2.67 

162.50 

24.41 

.66 

.20 

.25 

428. 

24.40 

1.33 

162.50 

18.32 

.95 

.40 

.71 

425. 

30.50 

1.07 

162.50 

15.29 

1.11 

.50 

1.04 

424. 

36.60 

.89 

162.50 

12.27 

1.34 

.60 

1.50 

424. 

42.70 

.76 

162.50 

9.31 

1.73 

.70 

2.27 

424. 

48.80 

.67 

162.50 

6.60 

2.53 

.80 

3.78 

426. 

54.90 

.60 

162.50 

5.80 

4.93 

.90 

8.29 

436. 

**»*♦»»**»  MINIMUM  COST  DESIGN  •♦••»»»» 

GATE 

FUNNEL 

DIMENSIONLESS 

TOTAL  COSTS 

PARAMETERS 

WIDTH  LENGTH 

VOLUME 

LENGTH 

PROJECTED 

b/B 

L/1 

LENGTH 

THOUSANDS  OF 

[M]  [M] 

[M^3] 

[M] 

[M] 

DOLLARS 

36.60  .89 

162.78 

12.27 

1.34 

.60 

1.50 

424. 

TOTAL  FUNNEL  WALL  AREA  [  M''2  ]  =  294.6 

TOTAL  GATE  WIDTH  WALL  AREA  [  M''2  ]  =  439.2 
TOTAL  GATE  LENGTH  WALL  AREA  [  M''2  ]  =  10.7 
TOTAL  FUNNEL/GATE  WALL  AREA  [  M''2  ]  =  744.4 
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Figure  9.  Funnel/Gate  Design  Example,  Selected  Contaminant  Fate  in  the  Gate. 
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B  =  61.m 
b  =  36.6  m 
{ -  .89  m 
L=  1.39  m 

The  last  four  lines  of  output  give  total  calculated  areas  for  the  fuimel  walls,  the  transverse  gate 
walls,  the  longitudinal  gate  walls,  and  the  summation  of  all  walls.  These  areas  can  be  combined 
as  needed  to  be  certain  estimated  system  areas  are  consistent  with  the  unit  costs  used  from  Table 
1 .  A  quick  view  of  Table  1  verifies  that  the  projected  area  of  744.4  m^  is  well  within  the 
specified  range  of  465  -  860  m^  used  to  select  values  for  unit  cost  coefficients,  C^avv,  Cjfg,,  and 
n 

'^sfgw* 
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SECTION  VI 


FGDM  CODE  LISTING 


Subroutine  FGDM.for: 


* 

*  FROM:  KIRK  HATFIELD 

*  DEPARTMENT  OF  CIVIL  ENGINEERING 

*  UNIVERSITY  OF  FLORIDA 

*  GAINESVILLE,  FLORIDA  32611  U.S.A. 

* 

************************************t***************^t************^****^:i::t:^.it:^ 

PROGRAM  FGDM 

* 

*  Variable  Declarations 

* 

c 

IMPLICIT  DOUBLE  PRECISION  (A-H,  0-Z) 

DIMENSION  BR(  1 0),GATEL(  1 0),VOL(  1 0),WING(1 0),RL(1 0), 

/TCOST(l  0),ANGLE(1 0),A(9,9),B(9,2),XX(4) 

COMMON  A,B,BR,GATEL, VOL, WING, RL,TCOST, ANGLE, XX 
READ  IN  THE  NAMES  OF  THE  INPUT  AND  OUTPUT  FILES 

CHARACTER*70INFIL,OUTFIL,INPUT,OUTPUT 
CHARACTER*!  HEADER 
WRITE(6,*)  'NAME  OF  INPUT  FILE* 

READ  (5,5)  INFIL 
INPUT=INFIL 

WRITE(6,*)  'NAME  OF  OUTPUT  FILE' 

5  FORMAT(A70) 

READ  (5,5)  OUTFIL 
OUTPUT  =  OUTFIL 

OPEN  (UNIT=1 1  ,FILE=INPUT,STATUS='OLD') 

OPEN  (UNIT=1 2,FILE=OUTPUT,STATUS='UNKNOWN') 
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o  n  o 


C  READING  IN  HEADER  (IGNORE) 

READ(11,70)HEADER 
70  FORMAT(Al) 

READ(11,70)HEADER 

C  READ  IN  WATER  QUALITY  STANDARDS  IN  (ug/1),  FOR  PCE  =  STDl,  TCE  =  STD2, 
C  DCE  =  STDS,  AND  VINYL  CHLORIDE  =  STD4 
C  READ  IN  TIME  STEP  INCREMENT  =  DT,  IN  DAYS 
C 

READ  (1 1,*)STD1,STD2,STD3,STD4,DT 
READ  (1 1,70)HEADER 
C 

C  THE  PERMEABILITY  OF  THE  AQUIFER  =  RKP,  IN  METERS/DAY 
C  SATURATED  THICKNESS  =  RPHIl,  IN  METERS 
C  THE  AQUIFER  GRADIENT  =  GRAD,  IN  METERS 
C  THE  PLUME  WIDTH  =  BIGB,  IN  METERS 
C 

READ(11,70)HEADER 
READ  (1 1,70)HEADER 
READ  (1 1,*)BIGB,GRAD,RPHI1,RKP 
READ(11,70)HEADER 
Q  =  -GRAD*BIGB*RKP*RPHI1 
C 

C  THE  POROSITY  OF  THE  GATE  =  THETA 
C  THE  RATIO  OF  AQUIFER  TO  GATE  PERMEABILITY  =  RATIOK 
C  THE  DEPTH  TO  AQUICLUDE  OR  DEPTH  OF  FUNNEL/GATE  WALLS 
C 

READ(11,70)HEADER 
READ  (1 1,70)HEADER 
READ  (1  l,*)THETA,RATIOK,DEPTH 
READ(11,70)HEADER 

THE  GATE  WIDTH  TO  PLUME  WIDTH  RATIOS,  b/B  EXAMINED  =  BR(I) 

READ(11,70)HEADER 
READ(11,70)HEADER 
READ(11,70)HEADER 
READ(11,70)HEADER 
READ  (1 1,*)BR(1),BR(2),BR(3),BR(4),BR(5) 

READ  (11,70)HEADER 
READ(11,70)HEADER 
READ  (1 1,70)HEADER 
READ  (1 1,*)BR(6),BR(7),BR(8),BR(9),BR(10) 

READ(11,70)HEADER 
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c 

C  THE  INITIAL  CONCENTRATION  OF  PCE  =  CIO  IN  ug/1,  AND  PCE  DEGRADATION 
C  RATE 

C  COEFFICIENTS  LAMBDA  la  =  RLMIA  AND  LAMBDA  lb  =  RLMIB,  IN  [  1/DAYS  ] 

C 

READ(11,70)HEADER 
READ(11,70)HEADER 
READ  (1 1,70)HEADER 
READ  (1 1,70)HEADER 
READ  (1 1,*)C10,RLM1A,RLM1B 
READ  (1 1,70)HEADER 
C 

C  THE  INITIAL  CONCENTRATION  OF  TCE  =  C20  IN  ug/1,  AND  TCE  DEGRADATION 
C  RATE 

C  COEFFICIENTS  LAMBDA  2a  =  RLM2A  AND  LAMBDA  2b  =  RLM2B,  IN  [  1/DAYS  ] 

C^ 

READ(11,70)HEADER 
READ(11,70)HEADER 
READ  (1 1,*)C20,RLM2A,RLM2B 
READ  (1 1,70)HEADER 
C 

C  THE  INITIAL  CONCENTRATION  OF  DCE  =  C30  IN  ug/1,  AND  DCE  DEGRADATION 
C  RATE 

C  COEFFICIENTS  LAMBDA  3a  =  RLM3A  AND  LAMBDA  3b  =  RLM3B,  IN  [  1/DAYS  ] 

O' 

READ  (1 1,70)HEADER 
READ  (1 1,70)HEADER 
READ  (1 1,*)C30,RLM3A,RLM3B 
READ(11,70)HEADER 
C 

C  THE  INITIAL  CONCENTRATION  OF  VINYL  CHLORIDE  =  C40  IN  ug/1,  AND 

C  THE  VINYL  CHLORIDE  DEGRADATION  RATE  COEFFICIENT  LAMBDA  4a  =  RLM3  A 
C  IN  [  1/DAYS  ]  ’ 

C 

READ(11,70)HEADER 
READ(11,70)HEADER 
READ  (1 1  ,*)C40,RLM4A 
READ(11,70)HEADER 
C 

C  READ  IN  COSTS  COEFFICIENTS: 

C  CSFFW  =  COST  PER  UNIT  AREA  OF  FUNNEL  WALL  MATERIAL 

C  CSFGL  =  COST  PER  UNIT  AREA  OF  LONGITUDINAL  GATE  WALL 

C  CSFGW  =  COST  PER  UNIT  AREA  OF  TRANSVERSE  GATE  WALL 
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o  o  o  n 


C  E 1 1  =  LAND  COSTS  PER  UNIT  AREA 

C  E12  =  REACTIVE  MEDIA  COSTS  PER  UNIT  VOLUME 

C 

READ(11,70)HEADER 
READ(11,70)HEADER 
READ(11,70)HEADER 
READ(1I,70)HEADER 
READ(11,70)HEADER 
READ(1L70)HEADER 
READ(11,70)HEADER 
READ(11,70)HEADER 
READ(11,70)HEADER 
READ  (11,*)  CSFFW,CSFGL,CSFGW,E1 1,E12 
C  E 1  =  COST  PER  UNIT  LENGTH  OF  FUNNEL  WING  WALL 

C  E2  =  COST  PER  UNIT  LENGTH  OF  GATE  LENGTH 

C  E3  =  COST  PER  UNIT  GATE  WIDTH 
El  =  CSFFW*DEPTH 
E2  =  CSFGL*DEPTH 
E3  =  CSFGW*DEPTH 
READ(11,70)HEADER 

READ  IN  INITIAL  ESTIMATES  OF  LAGRANGIAN  MULTIPLIERS  (USE  SUGGESTED 
VALUES  FROM  REPORT) 

READ  (1 1,70)HEADER 
READ  (1 1,70)HEADER 
READ(11,70)HEADER 
READ(11,70)HEADER 
READ  (11,*)  XX6,XX7,XX8,XX9 
C 

C 

C  FIRST  WE  SOLVE  THE  TRANSPORT  PROBLEM,  IN  WHICH  WE  DETERMINE  THE 
C  MINIMUM  HYDRAULIC  RESIDENCE 

C  TIME  NEEDED  TO  MEET  WATER  QUALITY  STANDARDS  AT  THE  GATE  EXIT 
C  THE  SUBROUTINE  RETURNS  WITH  A  VALUE  FOR  THE  MINIMUM  DETENTION 
C  TIME,  TD, 

C  NEEDED  TO  ACHIEVE  SUFFICIENT  DEGRADATION  SUCH  THAT  WATER 
C  QUALITY  STANDARDS 
C  ARE  MET 
C 

CALL  TRANS(STD  1  ,STD2,STD3 ,STD4,DT,C  1 0,RLM  1  A,RLM  1  B,C20,RLM2 A, 
/RLM2B,C30,RLM3A,RLM3B,C40,RLM4A,TD) 
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o  u 


c 

^  :i(  if:  9ic  3|c  :|t  ^  ^  sf: 

Q  ^c:ic9ic3ic^3|c:|c3H^9|c3icifc3ic3tcHc’»: 

C 

C  NOW  WE  SOLVE  FOR  VARIOUS  FUNNEL/GATE  DESIGN  CONFIGURATIONS  THAT 
C  PRODUCE  SUFFICIENT 

C  HYDRAULIC  RESIDENCE  TIME,  TD,  WITHIN  THE  REACTIVE  GATE.  IN  ADDITION, 
C  A  FIRST  ESTIMATE 

C  IS  FOUND  FOR  THE  LEAST  COST  DESIGN  IS  FOUND. 

C 

CALL  HYD(TD,BR,GATEL,VOL,WING,RL,RATIOK,RKP,RPHIl  ,Q,BIGB,TCOST 
/,THETA,E1,E2,E3,E1 1,E12,K) 

C 

Q  *H:**4i4ft********** 

C 

C  NOW  WE  DETERMINE  THE  FUNNEL/GATE  DESIGN  CONFIGURATION  OF 
C  MINIMUM  COST 
C 

CALL  OPTI(BIGB,TD,Q,THETA,RPHIl  ,RKP,RATIOK,BR(K),WING(K),GATEL(K), 
/RL(K),E1  ,E2,E3,E1 1  ,E12,XX6,XX7,XX8,XX9,DEPTH) 

91  PRINT*, THE  END' 

CLOSE(UNIT=ll) 

CLOSE(UNIT=12) 

STOP 

END 
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Subroutine  TRANS.for: 


C 

C  SUBROUTINE  TRANS  SIMULATES  TRANSPORT  THROUGH  THE  GATE  OF  A 
C  FUNNEL/GATE  SYSTEM 
C 

C  PARAMETERS  USED  BY  THE  SUBROUTINE  INCLUDE: 

C  WATER  QUALITY  STANDARDS  IN  (ug/I),  FOR  PCE  =  STDl ,  TCE  =  STD2, 

C  DCE  =  STDS,  AND  VINYL  CHLORIDE  =  STD4, 

C  TIME  STEP  INCREMENT  =  DX,  IN  DAYS 

C  THE  INITIAL  CONCENTRATION  OF  PCE  =  C 1 0,  IN  ug/1,  AND  PCE  DEGRADATION 
C  RATE 

C  COEFFICIENTS  LAMBDA  1  a  =  RLM 1 A  AND  LAMBDA  1  b  =  RLM 1 B,  IN  [  1  /DAY] 

C  THE  INITIAL  CONCENTRATION  OF  TCE  =  C20,  IN  ug/1  AND  TCE  DEGRADATION 
C  RATE 

C  COEFFICIENTS  LAMBDA  2a  =  RLM2A  AND  LAMBDA  2b  =  RLM2B,  IN  [1/DAY] 

C  THE  INITIAL  CONCENTRATION  OF  DCE  =  C30,  IN  ug/1  AND  DCE  DEGRADATION 
C  RATE 

C  COEFFICIENTS  LAMBDA  3a  =  RLM3 A  AND  LAMBDA  3b  =  RLM3B,  IN  [  1  /DAY] 

C  THE  INITIAL  CONCENTRATION  OF  VINYL  CHLORIDE  =  C40,  IN  ug/I, 

C  AND  THE  VINYL  CHLORIDE  DEGRADATION  RATE  COEFFICIENT  LAMBDA  4a  = 
C  RLM4A,  IN  [1/DAY] 

C 

C  THE  SUBROUTINE  RETURNS  WITH  A  VALUE  FOR  THE  CRITICAL  HYDRAULIC 
C  RESIDENCE  TIME  =  X,  IN  DAYS 

C  NEEDED  TO  ACHIEVE  SUFFICIENT  DEGRADATION  SUCH  THAT  WATER 
C  QUALITY  STANDARDS 
C  ARE  MET 

SUBROUTINE  TRANS(STD  1  ,STD2,STD3,STD4,DX,C  1 0,RLM  1  A,RLM  1 B, 
/C20,RLM2A,RLM2B,C30,RLM3A,RLM3B,C40,RLM4A,X) 

IMPLICIT  DOUBLE  PRECISION  (A-H,  O-Z) 

CHARACTER'S  ATM 
CHARACTER*8  ATM6 

CHARACTER*  1 3  ATM  1 , ATM2, ATMS , ATM4,ATM5 
CHARACTER*38  ATM7 
ATM  =TIME' 

ATM1='  Cl  ' 

ATM2='  C2  ' 

ATM3='  C3  ’ 

ATM4='  C4  ' 

ATM5='  [ug/1]  ’ 

ATM6='[  DAY  ]’ 

C  PARAMETER  RJ,RK,RG,RM,RF,RH,  AND  RP  ARE  SIMPLY  INTERMEDIATE 
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C  PARAMETERS  DEFINED 

C  IN  THE  THEORY  SECTION  OF  THE  REPORT  AS  j,k,g,m,f,h,p. 

RJ  =  (RLMIB+RLMI  A) 

RK  =  RLM1A 

RG  =  (RLM2A+RLM2B) 

RM  =  (RLM3A+RLM3B) 

RF  =  RLM2A 
RH  =  RLM4A 
RP  =  RLM3A 

C  CHECK  TO  MAKE  SURE  PARMETERS  ARE  NEVER  EQUAL 
C  OTHERWISE  A  MORE  EXPLICIT  SOUTION  MUST  BE  MADE  FOR  ALL 
C  POSSIBLE  COMBINATIONS  OF  EQUAL  PARAMETER  PAIRS  INCLUDING  RJ,RG,RM, 
C  AND  RH. 

RJ  =  RJ+.OOOOOl 
RG  =  RG+.000002 
RH  =  RH+.000003 
RM  =  RM+.000004 

C  SPECIFYING  THE  INITIAL  CONDITIONS  AT  THE  ENTRANCE  OF  THE  GATE  OR 
C  WHERE  THE  TOTAL  HYDRAULIC 

C  RESIDENCE  TIME,  X,  IS  ZERO,  WHILE  Cl,  C2,  C3,  AND  C4  RESPECTIVELY 
CORRESPOND  TO  DISSOLVED 

C  CONCENTRATION  (IN  ug/1)  VARIABLES  FOR  PCE,TCE,DCE,  AND  VINYL 
CHLORIDE 
X=0.0 
C1=C10 
C2=C20 
C3-C30 
C4=C40 

C  WRITE  OUT  THE  HEADING  FOR  THE  OUTPUT  GENERATED  BY  TRANS 
C 

WRITE(  1 2,2 1 )  ATM,  ATM  1  ,ATM2,  ATM3,ATM4 
WRITE(12,*)” 

WRITE(  1 2,20)ATM6,ATM5,ATM5,ATM5,ATM5 
WRITE(12,*)' ' 

WRITE(12,*)" 

20  F0RMAT(1X,A8,2X,4(13A)) 

21  FORMAT(A8,2X,4(13A)) 

C  WRITING  OUT  THE  INITIAL  CONDITIONS  AT  THE  EXIT  OF  THE  GATE 
C  WHICH  INCLUDE  THE  CUMULATIVE  HYDRAULIC  RESIDENCE  TIME,  X,  AND  THE 
C  FOUR  CONSTITUENT 

C  CONCENTRATIONS,  PCE,  TCE,  DCE,  AND  VINYL  CHLORIDE 
WRITE(12,68)  X,C1,C2,C3,C4 
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C  THIS  IS  THE  BEGINNING  OF  A  LOOP  THAT  SOLVES  FOR  THE  CONCENTRATION 
C  OF  EACH  DISSOLVED  CONSTITUENT  AS  THE  TOTAL  HYDRAULIC  RESIDENCE 
C  TIME,  X  IS  INCREASED  BY 

C  AN  INCREMENT,  DX,....THE  PROGRAM  CONTINUES  TO  SOLVE  FOR 
CONSTITUENT 
C  CONCENTRATIONS 

C  FOR  EVER  INCREASING  HYDRAULIC  RESIDENCE  TIMES,  X,  UNTIL  CALCULATED 
C  CONSTITUENT  CONCENTRATIONS 

C  ARE  LESS-THAN-OR-EQUAL-TO  PERTINENT  WATER  QUALITY  STANDARDS 
10  X  =  X+DX 

C  THE  SOLUTION  TO  THE  FIRST  REACTANT  PCE  (SEE  REPORT) 

Cl  =  C10*EXP(-RJ*X) 

C  THE  SOLUTION  TO  THE  SECOND  REACTANT  TCE  (SEE  REPORT) 

TERMl  =  RK*C10 

TERM2  =  (EXP(-RJ*X)-EXP(-RG*X))/(RG-RJ) 

C2  =  C20*EXP(-RG*X)+TERM1*TERM2 
C  THE  SOLUTION  TO  THE  THIRD  REACTANT  DCE  (SEE  REPORT) 

RAl  =  C30-RF*C20/(RM-RG) 

RA2  =  -RF*RK*C10/((RG-RJ)*(RM-RJ)) 

RA3  =  RF*RK*C10/((RG-RJ)*(RM-RG)) 

RA  =  RA1+RA2+RA3 
RDl  =  RF*C20/(RM-RG) 

RD2  =  -RF*RK*C10/((RG-RJ)*(RM-RG)) 

RD  =  RD1+RD2 

RE  =  RF*RK*C10/((RG-RJ)*(RM-RJ)) 
C3=RA*EXP(-RM*X)+RE*EXP(-RJ*X)+RD*EXP(-RG*X) 

C  THE  SOLUTION  TO  THE  FOURTH  REACTANT  VINYL  CHLORIDE  (SEE  REPORT) 
TERMl  =  C40*EXP(-RH*X) 

TERM2  =  RA*(EXP(-RM*X)-EXP(-RH*X))/(RH-RM) 

TERM3  =  RE*(EXP(-RJ*X)-EXP(-RH*X))/(RH-RJ) 

TERM4  =  RD*(EXP(-RG*X)-EXP(-RH*X))/(RH-RG) 

C4  =  TERMl +RP*(TERM2+TERM3+TERM4) 

C  WRITING  OUT  THE  CUMULATIVE  HYDRAULIC  RESIDENCE  TIME,  X,  AND  THE 
C  FOUR  CONSTITUENT 

C  CONCENTRATIONS,  PCE,  TCE,  DCE,  AND  VINYL  CHLORIDE 
WRITE(12,68)  X,C1,C2,C3,C4 
68  FORMAT(F6.2,4(3X,F10.2)) 

C  THIS  IS  WHERE  THE  PROGRAM  HAS  DETERMINES  IF  CALCULATED 
C  CONSTITUENT  CONCENTRATIONS 

C  ARE  LESS-THAN-OR-EQUAL-TO  PERTINENT  WATER  QUALITY  STANDARDS 
IF(Cl.GT.STDl)GOTO  10 
IF(C2.GT.STD2)  GO  TO  10 
IF(C3.GT.STD3)  GO  TO  10 
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o  o 


IF(C4.GT.STD4)  GO  TO  10 


C 

C 

C 

C 

C 

C 


THIS  IS  THE  END  OF  THE  LOOP  THAT  SOLVES  FOR  THE  CONCENTRATIONS 

™  hydraulic  residence 

AN  INCREMENT,  DX, 


WRITE(12,*)'  ' 

ATM7  ='  CRITICAL  HYDRAULIC  RESIDENCE  TIME  = ' 
WRITE(12,70)  ATM7,X,ATM6 
70  FORMAT(A38,2X,F6.2,2X,A8) 

WRJTE(12,*)" 


THE  PROGRAM  HAS  DETERMINED  THAT  CALCULATED 
CONCENTRATIONS 


CONSTITUENT  C  C 


c  m  pertinent  water  quality  standards  and 

C  MEET  S  S  ™ 


RETURN 

END 


54 


Subroutine  HYD.for: 


C 

C  SUBROUTINE  HYD  CONTAINS  APPROPRIATE  RELATIONSHIPS  TO  SOLVE  FOR 
C  VARIOUS  FUNNEL/GATE  DESIGN  CONFIGURATIONS  THAT  PRODUCE 
SUFFICIENT 

C  HYDRAULIC  RESIDENCE  TIME,  T,  WITHIN  THE  REACTIVE  GATE. 

C  IN  ADDITION,  THIS  SUBROUTINE  GENERATES  AN  INITIAL  ESTIMATE 
C  OF  THE  DESIGN  THAT  REPRESENT  THE  MINIMUM  COST  CONFIGURATION. 

C  THE  PERTINENT  SYSTEM  DIMENSION  OF  THE  MINIMUM  COST  CONFIGURATION 
C  ARE  STORE  IN  SEVERAL  COMMON  BLOCK  ARRAYS. 

C 

C  PARAMETERS  USED  BY  THE  SUBROUTINE  INCLUDE; 

C  THE  CRITICAL  HYDRAULIC  RESDENCE  TIME,  X,  IN  DAYS,  PROVIDED  BY 
C  SUBROUTINE  TRANS 

C  AN  ARRAY  CONTAINING  SEVERAL  GATE  WIDTH  TO  FUNNEL  WIDTH  RATIOS, 
C  b/B,  =  BR(I) 

C  THE  PERMEABILITY  OF  THE  AQUIFER=RKP,  IN  METERS/DAY 
C  SATURATED  THICKNESS=RPHI  1 ,  IN  METERS 

C  THE  AQUIFER  DISCHARGE=Q  (IN  CUBIC  METERS/DAY),  ACROSS  THE  FUNNEL 
C  WIDTH 

C  THE  FUNNEL  WIDTH=BIGB,IN  METERS 
C  THE  POROSITY  OF  THE  GATE  =  THETA 
C  THE  RATIO  OF  AQUIFER  TO  GATE  PERMEABILITY=RATIOK 
C  COSTS  COEFFICIENTS: 

C  El  =  COST  PER  UNIT  LENGTH  OF  FUNNEL  WING  WALL 

C  E2  =  COST  PER  UNIT  LENGTH  OF  GATE  LENGTH 

C  E3  =  COST  PER  UNIT  GATE  WIDTH 

C  El  1  =  LAND  COSTS  PER  UNIT  AREA 
C  E 1 2  =  COSTS  PER  UNIT  VOLUME  OF  REACTIVE  MEDIA 

C 

C  THE  SUBROUTINE  RETURNS  WITH  FUNNEL/GATE  DIMENSIONS  AND  COSTS 
C  STORED 

C  IN  SEVERAL  COMMON  BLOCK  ARRAYS  INCLUDING: 

C  GATEL(J)-AN  ARRAY  CONTAINING  CALCULATED  GATE  LENGTHS  (IN 
C  METERS) 

C  FOR  SYSTEM  CONFIGURATION  J 

C  VOL(J)-AN  ARRAY  CONTAINING  CALCULATED  GATE  VOLUMES  (IN  CUBIC 
C  METERS) 

C  FOR  SYSTEM  CONFIGURATION  J 

C  WING(J)-AN  ARRAY  CONTAINING  CALCULATED  FUNNEL  WALL  LENGTHS  (IN 
C  METERS) 

C  FOR  SYSTEM  CONFIGURATION  J 
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c  leSs  calculated  projected  funnel  wall 

C  FOR  SYSTEM  CONFIGURATION  J 

C  TCOST(J)~AN  ARRAY  CONTAINING  CALCULATED  TOTAL  SYSTEM  COSTS 
C  FOR  SYSTEM  CONFIGURATION  J 
C 


C  FINALLY  THE  SUBROUTINE  HYD  RETURNS  WITH  A  VALUE  FOR  K  AN  INTEGER 
C  VARIABLE  THAT  RECORDS  ’ 

C  WHICH  OF  THE  SEVERAL  SYSTEM  CONFIGURATIONS  EXAMINED  PRODUCED 
C  THE  LOWEST  CALCULATED  rKUL»uut,u 

c  mu^siom  ofTs  yste™  ^  ^ 

C  ™  Je  TOSroNS^®”"^'^  ™  CONFIGURATION. 

n  V  initial  decision  variable  estimate  in  subroutine  opti 

c  A  nonlinear  ’ 

CONHOU^TOn'^”^™^  designed  to  IDENTIFY  THE  TRUE  MINIMUM  COST 

c 


IMPLICIT  DOUBLE  PRECISION  (A-H,  0-Z) 

DIMENSION  BR(  1 0),GATEL(  1 0),  VOL(  1 0),  WING(  1 0),RL(  1 0) 
/TCOST(10),ANGLE(10)  v  y  v 

CHARACTER*6  ATM1,ATM2,ATM3 

CHARACTER*?  ATM4, ATMS, ATM7, ATMS, ATMl  1,ATM12  ATM19 
CHARACTER*9  ATM9 
CHARACTER*  13  ATM  10 
CHARACTER*  11  ATM13,ATM14 
CHARACTER*!  2  ATM6 
C 

C  FOR  VARIOUS  b/B  RATIOS  DETERMINE  THE  LENGTH  OF  THE 
C  REACTOR  ZONE  "1"  AND  THE  PROJECTED  FUNNEL  WALL  LENGTH  "L" 
C  (SEE  REPORT  FOR  DEFINITIONSC  OF  "L"  AND  "1"1 
C 

c 

C  OUTPUT  HEADER 
ATMl  = '  WIDTH' 

ATM2  =  'LENGTH' 


ATM3  =  'VOLUME' 

ATM4  = '[  M^3  ]' 

ATMS  = '  [  M  ] ' 

ATM6  =  'THOUSANDS  OF' 


S6 


n  n  n  n  n 


ATM7  = '  GATE  ' 

ATMS  = '  FUNNEL' 

ATM9  =  'PROJECTED' 

ATMIO  =  'DIMENSIONLESS' 

ATM  14  = '  PARAMETERS' 

ATM11='  b/B  ' 

ATM12  ='  L/1  ' 

ATM13  =  'TOTAL  COSTS' 

ATM19  =  'DOLLARS' 

C 

WRJTE(12,*)" " 

WRITE(12,72)  ATM7,ATM8,ATM10,ATM13 

72  FORMAT(10X,A7,16X,A7,1 1X,A13,5X,A1 1) 

WRITE(12,74)  ATM14 

74  F0RMAT(52X,A11) 

WRITE(12,73)  ATM1,ATM2,ATM3,ATM2,ATM9,ATM1  1,ATM12 

73  F0RMAT(1X,A6,3X,A6,3X,A6,4X,A6,3X,A9,2X,A7,3X,A7) 

WRITE(  12,75)  ATM2,ATM6 

75  FORMAT(40X,A6,23X,A12) 

WRITE(12,76)ATM5,ATM5,ATM4,ATM5,ATM5,ATM19 

76  F0RMAT(1X,A7,2X,A7,1X,A7,4X,A7,4X,A7,24X,A7) 

WRITE(12,*) 

C 

C  INITIAL  VALUE  OF  RLOLD  =  0.  AND  K  =  0  WHERE  K  IS  DEFINED  ABOVE  AND 
CRLOLD 

C  STORES  CALCULATED  VALUES  OF  THE  PROJECTED  FUNNEL  LENGTHS 
C  BETWEEN  SUCCESSIVE  INTERATIONS  USED  TO  OBTAIN  A  CONVERGENT 
C  SOLUTION  TO 

C  PROJECTED  FUNNEL  WALL  LENGTH  RL(J)  FOR  THE  SPECIFIC  b/B  RATIO  OR  BR(J) 
OF  INTEREST. 

C 

RLOLD  =  0. 

K=0 


THE  BEGINNING  OF  A  LOOP  THAT  CYCLES  THOUGH  10  POTENTIAL  VALUES  OF 
THE 

b/B  RATIO  TO  CALCULATED  10  POTENTIAL  FUNNEL/GATE  CONFIGURATIONS 

DO  90 1=1,10 
C 

C  TCOST  'TOTAL  COST  AN  INITIALLY  LARGE  VALUE  SO  THAT  THE  CHECK  IF 
C  STATEMENT  BEFORE 

C  THE  CONTINUE  STATEMENT  SEES  A  DECREASING  COST  AND  DOES  NOT 
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C  PREMATURELY  ASSIGN  K 
C  AN  INCORRECT  VALUE 
C 

TCOST(I)  =  1000000000000. 

c 

C  RATIOB  HOLD  THE  CURRENT  VALUE  OF  b/B  BEING  USED  TO  CONFIGURE  A 
C  SYSTEM 

RATIOB  =  BR(I) 

C 

C  FROM  CHRISTENSEN  AND  HATFIELD  (1994)  IT  IS  PREFERABLE  TO  NOT 
C  EXAMINE  SYSTEM  CONFIGURATIONS  WITH  b/B  RATIOS  LESS  THAN  OR  EQUAL 
C  TO  k/ko  =  RATIOK 

IF(BR(I).LE.RATIOK)  GO  TO  90 
C 

C  THE  FOLLOWING  CLUSTER  OF  5  LINES  GIVE  AN  INITIAL  CALCULATION  OF  THE 
C  PROJECTED 

C  LENGTH  OF  THE  FUNNEL  OR  WING  WALL  RL(I).  NOTE:  THE  DURING  THIS 
C  INITIAL  CALCULATION 

C  TERMS,  WHICH  IS  EQUAL  TO  THE  SATURATED  THICKNESS  OF  THE  AQUIFER  AT 
C  THE  GATE  ENTRANCE, 

C  IS  APPROXIMATED  WITH  RPHIl,  THE  SATURATED  AQUIFER  THICKNESS  AT  THE 
C  UPGRADIENT  ENTRANCE 
C  OF  THE  FUNNEL 
TERMS  =  RPHIl 

TERMl  =0.5*((RATIOK/RATIOB)-1.) 

TERM2  =  DLOG(RATIOB)/(l RATIOB)  +  1 . 

TERM4  =  Q*X/(THETA*BIGB*RATIOB) 

50  RL(I)  =  TERMl  *TERM4/(TERMS*TERM2) 

C 

C  AFTER  THE  FIRST  INTERATION,  TERMS  IS  NOW  CALCULATED  USING  A 
C  PREVIOUSLY  CALCULATED  RL(I) 

C  AS  DEFINED  BY  EQUATION  *****  IN  THE  REPORT 
C 

TERMS  =  (RPHI1*RPHI1  +  2.0*RL(I)*Q*DLOG(RATIOB)/ 

/(RKP*BIGB*(1 .0-RATIOB)))**0.5 
C 

C  ABSOLUTE  RELATIVE  DEFERENCE  BETWEEN  SUCCESSIVE  CALCULATIONS  OF 
C  RL(I)  ARE  EXPRESSED  WITH  THE 

C  VARIABLE  'ERROR'  WHERE  RLOLD  CONTAINES  THE  CALCULATED  VALUE  OF 
C  RL(I)  FROM  THE  PREVIOUS 
C  ITERATION 
C 

ERROR  =  ABS((RL(I)-RLOLD)/RL(I)) 
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RLOLD  =  RL(I) 

C  IF  THE  DIFFERENCE  BETWEEN  RL(I)  VALUES  CALCULATED  BETWEEN 
C  INTERATIONS  IS  GREATER  THAN  1  PERCENT 
C  THE  PROGRAM  RETURNS  TO  LINE  50  AND  RECALCULATES  RL(I) 
IF(ERROR.GT..0 1 )  GO  TO  50 
C 
C 

C  CALCULATED  A  GATE  LENGTH  =  GATEL(I) 

C 

GATEL(I)  =  RL(I)*TERM2/TERM1 
C 

C  CHECK  RESIDENCE  TIME,  X  IN  THE  REACTOR,  HERE  THE  RESIDENCE  TIME  IS 
C  RECALCULATED  TO  COMPARE 
C  WITH  THE  INITIALLY  SPECIFIED  VALUE 
C 

TERM4  =  Q/(THETA*BIGB*RATIOB) 

TERMS  =  (RPHI1*RPHI1  +  2.0*RL(I)*Q*DLOG(RATIOB)/ 

/(RKP*BIGB*(1 .0-RATIOB)))**0.5 
X  =  GATEL(I)*TERM3/TERM4 
C 

C  CALCULATE  GATE  VOLUME 
C 

VOL(I)  =  GATEL(I)*RATIOB*BIGB*TERM3 
C 

C  CALCULATE  LENGTH  OF  WING  WALL 
C 

TERM  =  (BIGB*BIGB*(1.-  RATIOB)**2.0)/4.0 
WING(I)  =  (RL(I)*RL(I)  +  TERM)**0.5 
C 

C  CALCULATE  ANGLE  OF  WING  WALL 
C 

ANGLE(I)  =  ASIN(BIGB*(l.-RATIOB)/(2.*WING(I))) 

C 

C  CALCULATE  TOTAL  FUNNEL/GATE  COST 
C 

TC0ST(I)=4.*E1*WING(I)+2.*(GATEL(I)*E2  +  E3*BIGB*RATIOB) 

/+E1 1  *BIGB*(GATEL(I)+2.*RL(I))+E12*VOL(I) 

C 

C  OUTPUT  STATMENTS  AS  NEEDED  TO  CREATE  A  TABLE  THAT  SHOWS  FOR 
C  SEVERAL  VALUES 

C  OF  THE  b/B  RATIO,  DIMENSIONS  OF  ASSOCIATED  FUNNEL/GATE  SYSTEM  AND 
C  THEIR  TOTAL  COST 
C 
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WRJTE(12,71)RATIOB*BIGB,GATEL(I),VOL(I),WING(I),RL(I),RATIOB 

/RL(I)/GATEL(I),TCOSTa) 

71  F0RMAT(1X,F6.2,2X,F6.2,4X,F6.2,2X,2X,F6.2,4X,F6.2,2X,F7.2,4X  F7  2 
/,2X,F10.0) 

C 

C  THE  FOLLOWING  CONDITIONAL  STATEMENT  DETERMINES  AT  WHICH 
C  ITERATION,  I,  A  DESIGN 

C  CONFIGURATION  WAS  FOUND  IN  WHICH  THE  TOTAL  COST  INCREASED  FROM 
C  THE  PREVIOUS 

C  ITERATION;  THUS,  DEFINING  THE  ITERATION  K  =  I-l,  WHERE  A  SPECIFIED  b/B 
C  RATIO  PRODUCED 

C  THE  LOWEST  COST  SYSTEM  CONFIGURATION. 

IF(LGT.  1  .AND.TCOST(I).GT.TCOST(I-1).AND.K.EQ.O)K=I-1 
90  CONTINUE 
C 

C  SUBROUTINE  HAS  COMPLETED  CALCULATIONS  FOR  A  SUITE  OF  SYSTEMS 
C  DEFINED  BY  THE  CRITICAL  HYDRAULIC 
C  RESIDENCE  TIME,X  AND  THE  b/B  RATIO 
RETURN 
END 
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Subroutine  OPTI.for: 


C  SUBROUTINE  OPTI  CONTAINS  APPROPRIATE  RELATIONSHIPS  THE  MINIMUM 
CCOST 

C  FUNNEL/GATE  DESIGN  CONFIGURATIONS  THAT  PRODUCE  SUFFICIENT 
C  HYDRAULIC  RESIDENCE  TIME,  TD,  WITHIN  THE  REACTIVE  GATE. 

C 

C  PARAMETERS  USED  BY  THE  SUBROUTINE  INCLUDE: 

C  THE  PLUME  WIDTH=BIGB,IN  METERS 

C  THE  MINIMUM  REQUIRED  HYDRAULIC  DETENTION  TIME,  X,  IN  DAYS, 

C  PROVIDED  BY 
C  SUBROUTINE  TRANS 

C  THE  AQUIFER  DISCHARGE=Q  (IN  CUBIC  METERS/DAY),  ACROSS  THE  PLUME 
C  WIDTH 

C  THE  POROSITY  OF  THE  GATE  =  THETA 

C  SATURATED  THICKNESS=RPHI  1 ,  IN  METERS 

C  THE  PERMEABILITY  OF  THE  AQUIFER=RKP,  IN  METERS/DAY 

C  THE  RATIO  OF  AQUIFER  TO  GATE  PERMEABILITY=RATIOK 

C  AN  INITIAL  ESTIMATE  OF  THE  GATE  WIDTH  TO  FUNNEL  WIDTH  RATIO,  b/B,  = 

CXI, 

C  GENERATED  BY  SUBROUTINE  HYD  AND  TRANSFERRED  THROUGH 

C  VARIABLE  XXI 

C  AN  INITIAL  ESTIMATE  OF  THE  FUNNEL  WALL  LENGTH  (IN  METERS),  =  X2, 

C  GENERATED  BY  SUBROUTINE  HYD  AND  TRANSFERRED  THROUGH 

C  VARIABLE  XX2 

C  AN  INITIAL  ESTIMATE  OF  THE  GATE  LENGTH  (IN  METERS),=  X3 , 

C  GENERATED  BY  SUBROUTINE  HYD  AND  TRANSFERRED  THROUGH 
C  VARIABLE  XX3 

C  AN  INITIAL  ESTIMATE  OF  THE  PROJECTED  FUNNEL  WALL  LENGTHS  (IN 
C  METERS),  =  X4, 

C  GENERATED  BY  SUBROUTINE  HYD  AND  TRANSFERRED  THROUGH 

C  VARIABLE  XX4 
C  COSTS  COEFFICIENTS: 

C  E 1  =  COST  PER  UNIT  LENGTH  OF  FUNNEL  WING  WALL 
C  E2  =  COST  PER  UNIT  LENGTH  OF  GATE  LENGTH 
C  E3  =  COST  PER  UNIT  GATE  WIDTH 
C  E 1 1  =  LAND  COSTS  PER  UNIT  AREA 

C  E 1 2  =  COSTS  PER  UNIT  VOLUME  OF  REACTIVE  MEDIA 

C  INITIAL  ESTIMATES  OF  LAGRANGIAN  MULTIPLIERS  X6,X7,X8,  AND  X9 
C  DEPTH  OF  FUNNEL/GATE  WALLS  OR  DEPTH  TO  AQUICLUDE,  DEPTH 
C 

C  THE  SUBROUTINE  FINDS  THE  DIMENSIONS  OF  THE  MINIMUM  COST 
C  FUNNEL/GATE  CONFIGURATION 
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n  n  n  n 


C  THEN  RECORDS  THOSE  DIMENSIONS  (i.e.,  GATE  LENGTH,  FUNNEL  WALL 
C  LENGTH, 

C  PROJECTED  FUNNEL  WALL  LENGTH,  GATE  WIDTH  TO  PLUME  WIDTH  RATIO, 
C  AND  THE  SATURATED 

C  AQUIFER  THICKNESS  IN  THE  GATE)  TO  THE  OUTPUT  FILE.  FINALLY  THE 
C  SUBROUTINE  RECORDS 
C  THE  TOTAL  COST  OF  THE  SYSTEM. 

C 

C 

SUBROUTINE  OPTI(BIGB,X,Q, THETA, RPHI 1  ,RKP,RATIOK,XX  1  ,XX2,XX3  ,XX4, 
/E1,E2,E3,E11,E12,X6,X7,X8,X9, DEPTH) 

IMPLICIT  DOUBLE  PRECISION  (A-H,  0-Z) 

DIMENSION  A(9,9),B(9,2),XX(4) 

COMMON  A,B 

CHARACTER*6  ATM1,ATM2,ATM3 

CHARACTER*?  ATM4,ATM5,ATM7,ATM8,ATM1 1,ATM12,ATM19 

CHARACTER*9  ATM9 

CHARACTER*  1 3  ATMl 0 

CHARACTER*  12  ATM6 

CHARACTER*  11  ATM13,ATM14 

CHARACTER*38  ATMl 5,ATM16,ATM1 7,ATM1 8 


OUTPUT  HEADER 
ATMl  =  ’  WIDTH' 

ATM2  =  'LENGTH' 

ATM3  =  'VOLUME' 

ATM4  = '[  M^3  ]' 

ATM5  = '  [  M  ] ' 

ATM6  =  'THOUSANDS  OF' 

ATM7  =  'GATE  ' 

ATM8  = '  FUNNEL' 

ATM9  =  'PROJECTED' 

ATM  10  =  'DIMENSIONLESS' 

ATM  14  =  '  PARAMETERS' 

ATM11='  b/B  ' 

ATMl 2-'  L/1  ' 

ATM  13  =  'TOTAL  COSTS' 

ATMl  5  =  'TOTAL  FUNNEL/GATE  WALL  AREA  [  M^2  ]  = ' 
ATMl 6  =  'TOTAL  FUNNEL  WALL  AREA  [  M"'2  ]  = ' 

ATMl 7  =  'TOTAL  GATE  LENGTH  WALL  AREA  [  M''^2  ]  = ' 
ATMl  8  =  'TOTAL  GATE  WIDTH  WALL  AREA  [  M^2  ]  =  ' 
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ATM19  =  'DOLLARS’ 


C 

WRITE(12,*)" " 

WRITE(12,*)" " 

WRITE(12,*)" " 

WRITE(12,*)"  **********  MINIMUM  COST  DESIGN  ********  " 

WRITE(12,*)’' " 

WRITE(12,*)" " 

WRITE(12,72)  ATM7,ATM8,ATM10,ATM13 

72  FORMAT(10X,A7,I6X,A7,1  1X,A13,5X,A1  1) 

WRITE(12,74)  ATM14 

74  F0RMAT(52X,A11) 

WRITE(12,73)  ATM1,ATM2,ATM3,ATM2,ATM9,ATM1  1,ATM12 

73  F0RMAT(IX,A6,3X,A6,3X,A6,4X,A6,3X,A9,2X,A7,3X,A7) 

WRITE(12,75)  ATM2,ATM6 

75  FORMAT(40X,A6,23X,A12) 

WRITE(  1 2,76)  ATMS,  ATM5,ATM4,  ATMS,  ATM5 ,  ATM  1 9 

76  F0RMAT(1X,A7,2X,A7,1X,A7,4X,A7,4X,A7,24X,A7) 

WRITE(12,*) 

C 

C 

C  PARAMETERS  NP,MP,N,AND  M  ARE  DIMENSIONS  USED  IN  CONSTRUCTING 
C  ARRAYS  USED 

C  IN  SUBROUTINE  GAUSSJ...A  GAUSS-JORDAN  SOLVER  PRESENTED  IN 
C  "NUMERICAL  RECIPIES  IN  FORTRAN"  BY  PRESS  et  al.  (1993) 

C 

C  ONE  NEED  ONLY  REMEMBER  THAT  NP=N  WHICH  IS  THE  NUMBER  OF 
C  VARIABLES  AND  THE 

C  NUMBER  OF  EQUATIONS  TO  BE  SOLVED  AS  GENERATED  FROM  FIRST 

CREATING 

C  THE  LAGRANGIAN 

C  AND  THEN  SUBSEQUENTLY  DEFINING  THE  NECESSARY  CONDITIONS  FOR  A 
C  STATIONARY  POINT 

C  FOR  THE  LAGRANGAJN....AND  M  IS  THE  NUMBER  OF  KNOWN  CONSTANTS 
C  DEFINING  THE  RIGHT-HAND-SIDE 
C  OF  EACH  EQUATION. 

C 

NP=9 

MP=2 

N=9 

M=1 

C 

C  PARAMETER  DEFINITIONS  AS  SPECIFIED  IN  THE  REPORT 
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non 


E4  =  THE  FUNNEL  WIDTH  "BIGB" 


E4  =  BIGB 


C 

C  E5  =  EQUIVALENT  TO  THE  KNOWN  PRODUCT  OF 

C  THICKNESS, 

C  RPHI,  AND  GATE  LENGTH,!. 

C 


GATE  WIDTH  'b'  FLOW 


E5=X*Q/THETA 

C 

E7=RPHI1**2. 

E8=2.*Q/(RKP*BIGB) 

E9  =  BIGB*BIGB/4. 

E10  =  RATIOK 
XX(1)  =  X9 
XX(2)=10.*X9 
XX(3)=  100.*X9 
XX(4)  =  500.  *X9 
DO  800  IK=1,4 

C  INITIAL  VARIABLE  ESTIMATES 
XI  =XX1 
X2  =  XX2 
X3=XX3 
X4  =  XX4 
X5  =  RPHIl 
X9=XX(IK) 

C  PRINT*  ,X1,X2,X3,X4,X5 
ITER  =  0 

C  *************************+******+*+******+,^**,^**,^,^,^*,^,^^^^^ 

C  SOLUTION^  application  OF  THE  NEWTON-RAPHSON 

C  ALGORITHM 
C 

C  DEFINING  DELTA  VALUES 
DXl  =  0. 

DX2  =  0. 

DX3  =  0. 

DX4=  0. 

DX5  =  0. 

DX6  =  0. 

DX7  =  0. 

DX8  =  0. 
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DX9  =  0. 

C  CALCULATE  UPDATED  PARAMETER  ESTIMATES  AND  FUNCTIONS  VALUES 
1  ITER=1+ITER 
XI  =X1  +DX1 
X2  =  X2  +  DX2 
X3  =  X3  +  DX3 
X4  =  X4  +  DX4 
X5  =  X5  +  DX5 
X6  =  X6  +  DX6 
X7  =  X7  +  DX7 
X8  =  X8  +  DX8 
X9  =  X9  +  DX9 

TMl  =  2.*E3*E4  +  X6*E4*X3*X5 
TM2  =  X7*(-X5*X5  +  E7  -  E8*X4/X1) 

TM3  =  -  X8*E9*(-2.  +  2.*X1) 

TM4  =  X9*(2.*X4*(DL0G(X1)  +  2.  -  2.*X1)) 

TM5  =  X9*(-X3*(-E10  -  1.  +  2.*X1)) 

FI  =  TM1+TM2+TM3+TM4+TM5 
F2  =  4.*E1  +X8*2.*X2 
TMl  =2.*E2  +  E11*E4  +  X6*E4*X1*X5 
TM2  =  -  X9*(E10*(1.-X1)  -  XI  +  X1*X1) 

F3  =  TM1+TM2 

TMl  =2.*E11*E4-2.*X8*X4 

TM2  =  X9*2.*(X1  *DLOG(Xl)  +  XI  -  X1*X1) 

F4  =  TM1+TM2 

F5  =  X6*E4*X1*X3  +  X7*2.*X5*(1.-  XI) 

F6  =  E4*X1*X3*X5-E5 

F7  =  X5*X5*(1.-X1)  -  E7*(1.-X1)  -  E8*X4*DL0G(X1) 

F8  =  X2*X2  -  E9*(l 2.*X1  +  XI *X1)  -  X4*X4 
TMl  =2.*X4’'(X1*DL0G(X1)  +  X1  -X1*X1) 

TM2  =  -  X3*(E10*(1.-X1)  -  XI  +  X1*X1) 

F9  =  TMl  +  TM2 

C  DEFINING  UPDATED  MATRIX  ELEMENTS  TO  BE  USED  IN  THE  GAUSS J  FOR 

C  SUBROUTINE. 

B(1,1)  =  -F1 
B(2,l)  =  -F2 
B(3,l)  =  -F3 
B(4,l)  =  -F4 
B(5,l)  =  -F5 
B(6,l)  =  -F6 
B(7,1)  =  .F7 
B(8,l)  =  -F8 
B(9,l)  =  -F9 
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C  402)^  derivatives  OF  F1,F2,F3,F4,F5,F6,F7,F8,AND  F9  (EQS.4-24  THROUGH 

C  TAKEN  WITH  RESPECT  TO  THE  NINE  VARIABLES 
C 


TMl  =X7*E8*X4/(X1*X1)-X8*E9*2. 

TM2  =  X9*(2.*X4/X1  -  4.*X4  -  2.*X3) 

A(1,1)  =  TM1+TM2 
A(l,2)  =  0. 

A(l,3)  =  X6*E4*X5  -  X9*(-E10  - 1.  +  2.*X1) 

A(l,4)  =  -X7*E8/XI  +  2.*X9*(DLOG(XI)  +  2.  -  2.*X1) 

A(l,5)  =  X6*E4*X3  -  2.*X7*X5 

A(1,6)  =  E4*X3*X5 

A(1,7)  =  -  X5*X5  +  E7  -  E8*X4/XI 

A(1,8)  =  -E9*(-2.+2.*X1) 

TMl  =  2.*X4*(DL0G(X1)  +  2.  -  2.*X1) 

TM2  =  -X3*(-E10- 1.  +2.*X1) 

A(l,9)  =  TMl  +  TM2 
A(2,l)  =  0. 

A(2,2)  =  2.*X8 
A(2,3)  =  0. 

A(2,4)  =  0. 

A(2,5)  =  0 
A(2,6)  =  0. 

A(2,7)  =  0. 

A(2,8)  =  2.*X2 
A(2,9)  =  0. 

A(3,I)  =  X6*E4*X5  -  X9*(-E10  -  1.  +  2.*X1) 

A(3,2)  =  0. 

A(3,3)  =  0. 

A(3,4)  =  0. 

A(3,5)  =  X6*E4*X1 
A(3,6)  =  E4*X1*X5 
A(3,7)  =  0. 

A(3,8)  =  0. 

A(3,9)  =  -E10*(1.-X1)  +  X1  -X1*X1 
A(4,l)  =  X9*2.*(DL0G(X1)  +  2.  -  2.*X1) 

A(4,2)  =  0. 

A(4,3)  =  0. 

A(4,4)  =  -2.*X8 
A(4,5)  =  0. 
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A(4,6)  =  0. 

A(4,7)  =  0. 

A(4,8)  =  -  2*X4 

A(4,9)  =  2.*(X1  *DL0G(X1)  +  XI  -  XI  *X1) 

A(5,l)  =  X6*E4*X3  -  X7*2.*X5 
A(5,2)  =  0. 

A(5,3)  =  X6*E4*X1 
A(5,4)  =  0. 

A(5,5)  =  2.*X7*(1.-X1) 

A(5,6)  =  E4*X1*X3 
A(5,7)  =  2.*X5*(1.-X1) 

A(5,8)  =  0. 

A(5,9)  =  0. 

A(6,1)  =  E4*X3*X5 
A(6,2)  =  0. 

A(6,3)  =  E4*X1*X5 
A(6,4)  =  0. 

A(6,5)  =  E4*X1*X3 
A(6,6)  =  0. 

A(6,7)  =  0. 

A(6,8)  =  0. 

A(6,9)  =  0. 

A(7,l)  =  -  X5*X5  +  E7  -  E8*X4/X1 
A(7,2)  =  0. 

A(7,3)  =  0. 

A(7,4)  =  -  E8*DL0G(X1) 

A(7,5)  =  2.*X5*(1.-X1) 

Ail, 6)  =  0. 

A(7,7)  =  0. 

A(7,8)  =  0. 

A(7,9)  =  0. 

A(8,l)  =  -E9*(-2.  +  2.*Xl) 

A(8,2)  =  2.*X2 
A(8,3)  =  0. 

A(8,4)  =  -2.*X4 
A(8,5)  =  0. 

A(8,6)  =  0. 

A(8,7)  =  0. 

A(8,8)  =  0. 

A(8,9)  =  0. 

A(9,l)  =  2.*X4*(DL0G(X1)  +  2.  -  2.*X1)  -  X3*(-E10  - 1.+  2.*X1) 
A(9,2)  =  0. 

A(9,3)  =-E10*(l.-Xl)  +  Xl  -X1*X1 
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A(9,4)  =  2.*(X1*DL0G(X1)  +  X1  -X1*X1) 

A(9,5)  =  0. 

A(9,6)  =  0. 

A(9,7)  =  0. 

A(9,8)  =  0. 

A(9,9)  =  0. 

C 

C  NOW  WE  CALL  GAUSSJ  BECAUSE  WE  WILL  BE 

C  SOLVING  FOR  THE  UPDATE  DELTA  VALUES  (WHICH  WILL  BRING  US  CLOSER  TO 
C  THE  LOWEST  COST  DESIGN 
C 

CALL  GAUSSJ(A,N,NP,B,M,MP) 

DXl  =B(1,1) 

DX2  =  B(2,1) 

DX3  =  B(3,1) 

DX4  =  B(4,1) 

DX5  =  B(5,1) 

DX6-B(6,1) 

DX7  =  B(7,1) 

DX8  =  B(8,1) 

DX9  =  B(9,1) 

C 

C  CHECK  TO  SEE  IF  VARIABLE  XI  IS  APPROACHING  A  VALUE  OF  1  ...IF  SO  WE 
C  CHOOSE  A  NEW 

C  INITIAL  VALUE  FOR  LAGRANGIAN  VARIABLE  X9 
C 

IF(X1.GT..95)  GO  TO  800 

C  CHECK  RELATIVE  CHANGES  IN  PARAMETER  ESTIMATES  BETWEEN 
C  INTERATIONS 

IF(DX1.LT..000000001)  GO  TO  300 
IF(ABS(DX1/X1).GE..005)  GO  TO  1 

300  IF(DX2.LT..000000001)  GO  TO  301 
IF(ABS(DX2/X2).GE..005)  GO  TO  1 

301  IF(DX3.LT..000000001)GOT0  302 
IF(ABS(DX3/X3).GE..005)  GO  TO  1 

302  IF(DX4.LT..000000001)GOT0  303 
IF(ABS(DX4/X4).GE..005)  GO  TO  1 

303  IF(DX5.LT..000000001)  GO  TO  304 
IF(ABS(DX5/X5).GE..005)  GO  TO  1 

304  IF(DX6.LT..000000001)GOT0  305 
IF(ABS(DX6/X6).GE..005)  GO  TO  1 

305  IF(DX7.LT..000000001)GOT0  306 
IF(ABS(DX7/X7).GE..005)  GO  TO  1 
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306  IF(DX8.LT..000000001)  GO  TO  307 
IF(ABS(DX8/X8).GE..005)  GO  TO  1 

307  IF(DX9.LT..000000001)GOT0  312 
IF(ABS(DX9/X9).GE..005)  GO  TO  1 
GO  TO  312 

C  PRINT  FINAL  OPTIMUM  PARAMETER  ESTIMATES 
800  CONTINUE 
C312  PRINT*, XI, X2,X3,X4,X5 
C  PRINT*  ,X6,X7,X8,X9,ITER 

312  TCOST  =  4.*E1*X2  +2.*E2*X3  +  2.*E3*E4*X1  +  El  1*E4*(X3+2.*X4) 
/+E12*E5 

C  TAREA  =  TOTAL  FUNNEL/GATE  WALL  AREA  [  M^2  ] 

C  TFWA  =  TOTAL  FUNNEL  WALL  AREA  [  M'^2  ] 

C  TGLA  =  TOTAL  GATE  LENGTH  WALL  AREA  [  M^2  ] 

C  TGWA  =  TOTAL  GATE  WIDTH  WALL  AREA  [  M^2  ] 

TAREA  =  4.*DEPTH*X2  +2.*DEPTH*X3  +  2.*DEPTH*BIGB*X1 
TFWA  =  4.*DEPTH*X2 
TGLA  =  2.*DEPTH*X3 
TGWA  =  2.*DEPTH*BIGB*X1 

WRITE(6,*)"  **********  MINIMUM  COST  DESIGN  ********  " 
WRITE(6,*)" " 

WRITE(6,*)"  " 

WRITE(6,72)  ATM7,ATM8,ATM10,ATM13 
WRITE(6,74)  ATM14 

WRITE(6,73)  ATM1,ATM2,ATM3,ATM2,ATM9,ATM1 1,ATM12 
WRITE(6,75)  ATM2,ATM6 

WRITE(6,76)ATM5,ATM5,ATM4,ATM5,ATM5,ATM  1 9 
WRITE(6,*)"  " 

WRITE(6,71)X1*BIGB,X3,X1*BIGB*X5*X3,X2,X4,X1,X4/X3,TC0ST 
WRITE(6,*)"  " 

WRITE(6,*)"  " 

WRITE(6,*)"  " 

WRITE(6,78)ATM  1 5,TAREA 
WRITE(6,*)"  " 

WRITE(6,*)"  " 

WRITE(6,*)"  " 

78  FORMAT(2X,A38,F6.1) 

WRITE(12,71)X1*BIGB,X3,X1*BIGB*X5*X3,X2,X4,X1,X4/X3,TC0ST 
71  FORMAT(lX,F6.2,2X,F6.2,4X,F6.2,2X,2X,F6.2,4X,F6.2,2X,F7.2,4X,F7.2 
/,2X,F10.0) 

WRITE(12,*)"  " 

WRITE(12,*)"  " 

WRITE(12,*)"  " 
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WRITE(12,78)ATM16,TFWA 

WRITE(12,78)ATM1 8,TGWA 

WRITE(  1 2,78)ATM  1 7,TGL  A 

WRITE(12,78)ATM15,TAREA 

RETURN 

END 


Subroutine  GAUSSJ.for: 


SUBROUTINE  GAUSSJ(A,N,NP,B,M,MP) 

PARAMETER  (NMAX=50) 

IMPLICIT  DOUBLE  PRECISION  (A-H,  0-Z) 

DIMENSION  A(NP,NP)3(NP,MP),IPIV  (NMAX),INDXR(NMAX),INDXC(NMAX) 
DOll  J=1,N 
IPIV(J)=0 

1  CONTINUE 
DO  22 1=1, N 

BIG=0. 

DO  13  J=1,N 
IF(IPIV(J).NE.1)THEN 
D0  12K=1,N 
IF(IPIV(K).EQ.O)  THEN 
IF(ABS(A(J,K)).GE.BIG)  THEN 
BIG=ABS(A(J,K)) 

IROW=J 

ICOL=K 

ENDIF 

ELSE  IF  (IPIV(K).GT.l)  THEN 
PAUSE  'SINGULAR  MATRIX' 

ENDIF 

2  CONTINUE 
ENDIF 

3  CONTINUE 
IPIV(ICOL)=IPIV(ICOL)+l 
IF(IROW.NE.ICOL)  THEN 

DO  14L=1,N 
DUM=A(IROW,L) 

A(IROW,L)=A(ICOL,L) 

A(ICOL,L)=DUM 

4  CONTINUE 
D0  15L=1,M 

DUM=B(IROW,L) 

B(IROW,L)=B(ICOL,L) 

B(ICOL,L)=DUM 

5  CONTINUE 
ENDIF 

INDXR(I)=IROW 

INDXC(I)=ICOL 

IF  (A(ICOL,ICOL).EQ.O.)  PAUSE  'SINGULAR  MARTRIX' 

PIVINV=1  ./A(ICOL,ICOL) 
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A(IC0L,IC0L)=1. 

DO  16L=1,N 

A(ICOL,L)=A(ICOL,L)*PIVINV 

16  CONTINUE 
DO  17L=1,M 

B(ICOL,L)=B(ICOL,L)*PIVINV 

17  CONTINUE 
DO  21  LL=1,N 

IF(LL.NE.ICOL)THEN 

DUM=A(LL,ICOL) 

A(LL,ICOL)=0. 

DO  18L=1,N 

A(LL,L)=A(LL,L)-A(ICOL,L)*DUM 

18  CONTINUE 
DO  19L=1,M 

B(LL,L)=B(LL,L)-B(ICOL,L)*DUM 

19  CONTINUE 
ENDIF 

21  CONTINUE 

22  CONTINUE 
DO  24L=N,1,-1 

IF(INDXR(L).NE.INDXC(L))  THEN 
DO  23  K=1,N 
DUM=A(K,INDXR(L)) 
A(K,INDXR(L))=A(K,INDXC(L)) 
A(K,INDXC(L))=DUM 

23  CONTINUE 
ENDIF 

24  CONTINUE 
RETURN 
END 
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Example  Input  File: 


STDl  [ug/1]  STD2  [ug/1]  STD3  [ug/1]  STD4  [ug/1]  DT(TIMESTEP)  [DAYS] 

♦ 

5.0  5.0  70.0  1.0  .1 

* 

FUNNEL  WIDTH  [M]  GRADIENT  AQU.  THICKNESS  [M]  AQU.  HYDRAULIC  COND. 

[M/d] 

* 

61.  -0.005  5.  4.32 

* 

GATE  POROSITY  HYDRAULIC  COND.  RATIO  (AQUIFER/GATE)  DEPTH  OF 

SYSTEM  WALLS  [M] 

* 

.45  .1  6. 

* 

VARIOUS  RATIOS  OF  GATE  WIDTH  TO  PLUME  WIDTH  b/B  EXAMINED 
* 

BR(1)  BR(2)  BR(3)  BR(4)  BR(5) 

* 

.02  .05  .1  .2  .4 

* 

BR(6)  BR(7)  BR(8)  BR(9)  BR(IO) 

* 

.5  .6  .7  .8  .9 

* 

****DATA  PLUME  CONTAMINANT  CONCENTRATIONS  AND  DEGRADATION 

PARAMETERS 
♦ 

PLUME  Cl  CONC.  [ug/1] 

* 

5000.  36.9 

* 

PLUME  C2  CONC.  [ug/1] 

♦ 

5000.  3.69 

♦ 

PLUME  C3  CONC.  [ug/1] 

* 

0.  .8318 

* 

PLUME  C4  CONC.  [ug/1] 


LAMBDA  lA  [1/DAY]  LAMBDA  IB  [1/DAY] 
0. 

LAMBDA  2A  [1/DAY]  LAMBDA  2B  [1/DAY] 
33.18 

LAMBDA  3A  [1/DAY]  LAMBDA  3B  [1/DAY] 
0. 

LAMBDA  4A  [1/DAY] 
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0. 


.8318 


* 


*  READ  IN  COSTS  COEFFICIENTS:  [  IN  UNITS  OF  THOUSANDS  OF  DOLLARS  ] 

*  CSFFW  =  COST  PER  UNIT  AREA  OF  FUNNEL  WALL  MATERIAL 

*  CSFGL  =  COST  PER  UNIT  AREA  OF  LONGITUDINAL  GATE  WALL 

*  CSFGW  =  COST  PER  UNIT  AREA  OF  TRANSVERSE  GATE  WALL 

*  E 1 1  =  LAND  COSTS  PER  UNIT  AREA 

*  E 1 2  =  REACTIVE  MEDIA  COSTS  PER  UNIT  VOLUME 

* 

CSFFW  CSFGL  CSFGW  Ell  E12 

* 

.193  .193  .193  .001  1.725 

* 

*  READ  IN  INITIAL  ESTIMATES  OF  LAGRANGIAN  MULTIPLIERS,  X6,  X7,  X8,  AND  X9 

X6  X7  X8  X9 

♦ 

1.  1.  1.  10. 


74 


Example  Output  File: 


TIME 

Cl 

C2 

C3 

C4 

[DAY] 

[ug/1] 

[ug/1] 

[ug/1] 

[ug/1] 

.00 

5000.00 

5000.00  .00 

.00 

.10 

124.86 

586.66 

879.91 

OO 

.20 

3.12 

26.22 

873.98 

115.50 

.30 

.08 

.94 

806.89 

173.34 

.40 

.00 

.03 

742.58 

221.27 

.50 

.00 

.00 

683.31 

260.45 

.60 

.00 

.00 

628.77 

291.96 

.70 

.00 

.00 

578.59 

316.79 

.80 

.00 

.00 

532.41 

335.79 

.90 

.00 

.00 

489.91 

349.74 

1.00 

.00 

.00 

450.81 

359.32 

1.10 

.00 

.00 

414.83 

365.15 

1.20 

.00 

.00 

381.72 

367.76 

1.30 

.00 

.00 

351.25 

367.62 

1.40 

.00 

.00 

323.22 

365.16 

1.50 

.00 

.00 

297.42 

360.76 

1.60 

.00 

.00 

273.68 

354.73 

1.70 

.00 

.00 

251.84 

347.36 

1.80 

.00 

.00 

231.74 

338.92 

1.90 

.00 

.00 

213.24 

329.60 

2.00 

.00 

.00 

196.22 

319.62 

2.10 

.00 

.00 

180.56 

309.13 

2.20 

.00 

.00 

166.15 

298.27 

2.30 

.00 

.00 

152.89 

287.18 

2.40 

.00 

.00 

140.68 

275.96 

2.50 

.00 

.00 

129.46 

264.71 

2.60 

.00 

.00 

119.12 

253.49 

2.70 

.00 

.00 

109.62 

242.37 

2.80 

.00 

.00 

100.87 

231.42 

2.90 

.00 

.00 

92.82 

220.67 

3.00 

.00 

.00 

85.41 

210.16 

3.10 

.00 

.00 

78.59 

199.92 

3.20 

.00 

.00 

72.32 

189.98 

3.30 

.00 

.00 

66.55 

180.35 

3.40 

.00 

.00 

61.23 

171.05 

3.50 

.00 

.00 

56.35 

162.09 
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3.60 

3.70 

3.80 

3.90 
4.00 

4.10 

4.20 

4.30 

4.40 

4.50 

4.60 

4.70 

4.80 

4.90 
5.00 

5.10 

5.20 

5.30 

5.40 

5.50 

5.60 

5.70 

5.80 

5.90 
6.00 

6.10 

6.20 

6.30 

6.40 

6.50 

6.60 

6.70 

6.80 

6.90 
7.00 
7.10 
7.20 

7.30 

7.40 

7.50 
7.60 

7.70 
7.80 


.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 


.00 

51.85 

153.46 

.00 

47.71 

145.18 

.00 

43.90 

137.25 

.00 

40.40 

129.65 

.00 

37.17 

122.40 

.00 

34.21 

115.47 

.00 

31.48 

108.88 

.00 

28.97 

102.59 

.00 

26.65 

96.62 

.00 

24.53 

90.95 

.00 

22.57 

85.57 

.00 

20.77 

80.47 

.00 

19.11 

75.63 

.00 

17.58 

71.06 

.00 

16.18 

66.73 

.00 

14.89 

62.65 

.00 

13.70 

58.79 

.00 

12.61 

55.14 

.00 

11.60 

51.71 

.00 

10.68 

48.47 

.00 

9.82 

45.42 

.00 

9.04 

42.54 

.00 

8.32 

39.84 

.00 

7.65 

37.30 

.00 

7.04 

34.91 

.00 

6.48 

32.66 

.00 

5.96 

30.55 

.00 

5.49 

28.57 

.00 

5.05 

26.71 

.00 

4.65 

24.96 

.00 

4.28 

23.32 

.00 

3.93 

21.79 

.00 

3.62 

20.35 

.00 

3.33 

19.00 

.00 

3.07 

17.74 

.00 

2.82 

16.56 

.00 

2.60 

15.46 

.00 

2.39 

14.42 

.00 

2.20 

13.45 

.00 

2.02 

12.55 

.00 

1.86 

11.70 

.00 

1.71 

10.91 

.00 

1.58 

10.17 
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7.90 

.00 

.00 

1.45 

9.48 

8.00 

.00 

.00 

1.33 

8.83 

8.10 

.00 

.00 

1.23 

8.23 

8.20 

.00 

.00 

1.13 

7.67 

8.30 

.00 

.00 

1.04 

7.14 

8.40 

.00 

.00 

.96 

6.65 

8.50 

.00 

.00 

.88 

6.19 

8.60 

.00 

.00 

.81 

5.77 

8.70 

.00 

.00 

.75 

5.37 

8.80 

.00 

.00 

.69 

5.00 

8.90 

.00 

.00 

.63 

4.65 

9.00 

.00 

.00 

.58 

4.33 

9.10 

.00 

.00 

.53 

4.03 

9.20 

.00 

.00 

.49 

3.75 

9.30 

.00 

.00 

.45 

3.48 

9.40 

.00 

.00 

.42 

3.24 

9.50 

.00 

.00 

.38 

3.01 

9.60 

.00 

.00 

.35 

2.80 

9.70 

.00 

.00 

.32 

2.61 

9.80 

.00 

.00 

.30 

2.42 

9.90 

.00 

.00 

.27 

2.25 

10.00 

.00 

.00 

.25 

2.09 

10.10 

.00 

.00 

.23 

1.95 

10.20 

.00 

.00 

.21 

1.81 

10.30 

.00 

.00 

.20 

1.68 

10.40 

.00 

.00 

.18 

1.56 

10.50 

.00 

.00 

.17 

1.45 

10.60 

.00 

.00 

.15 

1.35 

10.70 

.00 

.00 

.14 

1.25 

10.80 

.00 

.00 

.13 

1.16 

10.90 

.00 

.00 

.12 

1.08 

11.00 

.00 

.00 

.11 

1.00 

11.10 

.00 

.00 

.10 

.93 

CRITICAL  HYDRAULIC  RESIDENCE  TIME  =  1 1 . 10  [  DAY  ] 
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GATE 

WIDTH  LENGTH  VOLUME 
[  M  ]  [  M  ]  [  M^3  ] 


FUNNEL  DIMENSIONLESS 
PARAMETERS 

LENGTH  PROJECTED  b/B  L/I 
LENGTH 
[M]  [M] 


total  costs 


THOUSANDS  OF 
DOLLARS 


12.20 

2.67 

162.50 

24.41 

24.40 

1.33 

162.50 

18.32 

30.50 

1.07 

162.50 

15.29 

36.60 

.89 

162.50 

12.27 

42.70 

.76 

162.50 

9.31 

48.80 

.67 

162.50 

6.60 

54.90 

.60 

162.50 

5.80 

.66 

.20 

.25 

428. 

.95 

.40 

.71 

425. 

1.11 

.50 

1.04 

424. 

1.34 

.60 

1.50 

424. 

1.73 

.70 

2.27 

424. 

2.53 

.80 

3.78 

426. 

4.93 

.90 

8.29 

436. 

**********  MINIMUM  COST  DESIGN  ******** 


GATE 


FUNNEL 


WIDTH 

LENGTH 

VOLUME 

LENGTH 

[M] 

[M] 

[M^3] 

[M] 

36.60 

.89 

162.78 

12.27 

DIMENSIONLESS  TOTAL  COSTS 
PARAMETERS 


LENGTH 

[M] 

1.34 


.60  1.50 


THOUSANDS  OF 
DOLLARS 


424. 


TOTAL  FUNNEL  WALL  AREA  [  M'^2  ]  =  294  6 

TOTAL  GATE  WIDTH  WALL  AREA  [  M^2  ]  =  439  2 
TOTAL  GATE  LENGTH  WALL  AREA  [  M^2  ]  =  107 
TOTAL  FUNNEL/GATE  WALL  AREA  [  M^2  ]  =  744  4 
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SECTION  VII 


CONCLUSIONS 

A  series  of  equations  which  constitute  an  analytical,  multi-contaminant,  degradation  and 
transport  model  have  been  developed.  This  model  can  be  applied  in  iterative  manner  to 
determine  the  critical  hydraulic  residence  time,  T,  required  to  achieve  specified  concentrations 
of  contaminants  within  a  gate  (e.g.,  to  be  certain  that  the  water  quality  at  the  gate  exit  meets  a 
specified  water  quality  standard  i.e.,  1  pg/l  vinyl  chloride).  To  use  this  transport  model,  the 
following  data  are  needed: 

a)  degradation  rate  coefficients  for  all  reactions  depicted  in  Figure  3; 

b)  specified  water  quality  goals  or  system  performance  standards  for  C,,  C2,  C3,  and 
C4  (e.g.,  primary  drinking  water  standards  for  PCE,  TCE,  DCE,  and  vinyl 
chloride);  and 

c)  the  concentration  of  each  contaminant  at  the  entrance  of  the  gate,  Cjo,  C20,  C30, 

and  C40.  Once  the  critical  hydraulic  residence  time  is  known,  the  hydraulic  model 
can  be  used  to  dimension  the  system. 

From  an  existing  analytical  hydraulic  model  for  fiiimel/gate  systems  developed  by 
Christensen  and  Hatfield  (1994),  an  iterative  design  algorithm  was  developed  for  dimensioning  a 
funnel/gate  system.  The  design  algorithm  requires  the  following  site  and  system  data  to  execute: 

a)  the  critical  hydraulic  residence  time  in  the  gate; 

b)  the  transverse  width  of  groundwater  flow  to  be  intercepted  within  the  funnel,  B; 

c)  the  saturated  thickness  of  the  phreatic  aquifer,  (|),  at  the  entrance  of  the  funnel; 

d)  the  groundwater  flow  through  the  area  defined  by  width  B  and  flow  thickness,  (J)]; 

e)  the  aquifer  hydraulic  conductivity; 

f)  the  porosity  of  the  reactive  porous  media  within  the  gate;  and 

g)  the  hydraulic  conductivity  of  the  porous  media  within  the  gate. 

Results  of  two  FRAC-3D  numerical  validation  studies  were  also  presented.  For  the  two 
funnel/gate  systems  examined,  the  numerical  results  showed  the  analytical  design  model 
produced  system  configurations  which  exhibited  80  to  85  percent  capture  efficiency.  These 
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efficiencies  increase  as  the  gate  becomes  wider.  Thus,  results  of  the  numerical  validation  suggest 
the  furmel/gate  design  model  could  be  used  to  pre-dimension  a  system;  however,  subsequent 
numerical  simulations  should  be  conducted  to  verify  the  hydraulics  of  the  design. 

A  funnel/gate  cost  estimation  model  was  also  developed  to  identify  minimum  cost 
designs  subject  to  four  constraints:  one  that  defines  the  desired  hydraulic  retention  time  in  the 
gate,  one  that  defines  the  funnel  wall  length  in  terms  of  other  system  dimensions,  and  two  that 
serve  as  hydraulic  constraints.  Costs  considered  in  the  model  include  land  costs,  costs  per  unit 
length  of  funnel  walls,  costs  per  unit  length  of  gate,  costs  per  unit  width  of  gate,  and  costs 
associated  with  the  reactive  medium  used  in  the  gate. 

The  general  cost  model  was  recast  into  a  Lagrangian  optimization  model.  As  a  result  of 
this  reformulation,  the  optimization  problem  was  reduced  to  that  of  solving  a  system  of  nine 
nonlinear  equations  with  nine  unknowns.  Recommendations  were  given  as  to  how  the  nine 
equations  could  be  solved. 

Finally  a  computer  program  was  developed  to  demonstrate  the  theory  that  was  developed. 
The  Funnel/Gate  Design  Model  (FGDM)  software  is  a  FORTRAN  program  developed  to 
ascertain  feasible  dimensions  of  funnel/gate  systems  and  also  identify  the  lowest  cost  design. 
Feasible  is  defined  here  to  include  designs  that  are  consistent  with  the  transport  and  hydraulic 
theories  outlined  in  Sections  2  and  3.  Examples  of  input  data  and  computational  results  are 
presented. 
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